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I. Progress Summary 


Work performed during the past year is summarized in this section. The discussion is subdivided into 
subject areas, as indicated. 

The following work has been completed in the area of solar and stellar magnetic fields, related photos- 
pheric phenomena and the relationships between magnetism, rotation, coronal and chromospheric 
emission in solar-like stars: 

A paper on the comparison between x-ray bright points (XBP) and He I 10830 dark points 
(Golub, Harvey, Herant and Webb 1989) was completed. Using several sets of simultaneous x-ray and 
He I 10830 images we compared dark points in He I with XBP on the x-ray images. For the largest 
and most obvious features there is a strong correlation; however, about two-thirds of the x-ray features 
were not identified on the He I images. Once an x-ray feature is located, it is nearly always possible to 
find a corresponding He I feature. The conclusion is that, using current selection methods, the He I 
10830 data cannot be used as a one-to-one proxy for the x-ray data. 

Schrijver 2et al. (1989) have derived an improved relation between magnetic flux (<S> B and chro- 
mospheric Ca II H+K line core emission flux densities (AF^j) for solar active and quiet regions. 
They find that approximates in contrast with previous determinations of a linear rela- 

tionship (which were biased due to analysis of data with a much smaller range in AF^j amd V 
The new AF^g-Og relation is nearly identical with a preliminary one found for less active dwar stars 
(Saar and Schrijver 1987), which in turn is consistent with stellar X-ray - magnetic flux and stellar 
flux-flux (AF^aii'Fx^ correlations (Saar and Schrijver 1987; Schrijver 1983). The agreement between 
relations derived for individual solar active regions and from hemisphere-averaged stellar data implies 
that important physical similarities exist between magnetic structures on the Sun and magnetically 
"quiet" dwaraf stars, at least at the chromospheric level. The internal consistency of the solar and stel- 
lar chromospheric and coronal relations further implies that the same may hold true in stellar coronae. 
The nonlinearity of the stellar AF ^^ -Fy may then suggest that the chromospheres and coronae of 
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cool stars are heated by different mechanisms. 

Bopp et al. (1989) have made a detailed study of the interesting binary star, HD 17433. They 
find that the star lies above the main sequence and may be quite young (approximately Pleiades age). 
Magnetic activity levels in HD 17433 are high, with chromospheric, transition region, and coronal 
emission fluxes at the stellar suraface up to 200 times that of the active Sun. Consistent with these 
elevated activity levels, Bopp et al. (using an improved magnetic line profile model) estimate that 
roughly 66% of the photosphere is covered by magnetic regions with a mean field strength of 2000 
Gauss. This represents the irst successful detection of a magnetic field on an active subgiant. Interest- 
ingly, the derived magnetic parameters are considerably larger than those expected from extrapolation 
of the preliminary relations found for dwarafs (Saar 1987). 

The TiO bands are sensitive indicators of temperature in stellar atmospheres. The very presence 
of these features in dwarf stars with T eff > 4250 K (approx K5V) is direct evidence for cool sp[ots. 
Saar and Nefff (1988) used observations of two TiO bands with differing temperature sensitivities to 
extract both the temperature (from the relative band depths) and the projected area (from the absolute 
band strengths relative to the continuum) of stellar spots. Two X-ray bright K dwarfs were studied 
and spots with AT approx 750 K were found to cover approximately 20% of the surface of BD+26730. 
The spectroscopic technique of determining star spot properties is superior to photometric techniques, 
which have difficulty detecting spots if there is no rotational modulation o the stellar light. A full 
description of this work will be submitted to ApJ. Letters shortly. 

Theoretical studies have been carried out in the following areas, among others : 

1. Neutral Beams as the Dominant Energy Transport Mechanism in Two-Ribbon Flares 

We have presented an alternative to die model of mildly relativistic electron beam for the tran- 
sportation of flare-energy from the corona to the chromosphere (see Dennis (Sol. Phys. 100): neutral 
berams of protons and electrons with equal velocity and with the protons, by virtue of their larger 
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mass, carrying most of the beam energy. In Martens (1988) the basic physical model is derived and 
found to be in agreement with observational requirements. In Martens and Young (1989) we make a 
detailed comparison with the observationally much better documented acceleration of protons in the 
geomagnetic tail, and find great similarity between the two. In Korevaar and Martens (1989) recent 
observations of polarized emission lines are reviewed, which clearly favor neutral beams over electron 
beams. Also, the dimensionless model equations for the proton acceleration are derived and found to 
have the same parameters as those for the geotail: this makes numerical results for the geotail directly 
applicable to solar two-ribbon flares. 

2. MHD and Circuit Models for Filament Eruptions 

We have modelled the formation and eruption of solar prominences with different approaches and 
compared the results. An MHD equilibrium model (van Ballegooijen and Martens 1989) has demon- 
strated the formation of a filament as the result of cancellation of photospheric flux at the neutral line, 
and the subsequent loss of equilibrium as more and more flux cancels. This is consistent with the 
scenario derived from recent observations by Martin and coworkers. A simplified circuit code (Mar- 
tens and Kuin 1989) reproduces the results of the MHD equilibrium calculations exactly and moreover, 
is also able to follow the filament after its loss of equilibrium. The filament is found to erupt, and its 
velocity curve is in agreement with that of typical filament eruptions. We are now modeling specific 
flares for which detailed velocity curves are available from observations. 

3. Studies of Radio Emission Mechanisms in Transient Events 

Tearing modes and Alfven waves are important phenomena in plasma astrophysics and in labora- 
tory plasma. The periodicity of tearing modes has been observed in experiments of tokamak plasma 
(Sunnji Tsiji, 1983). Huang (1987) derived a complex eigenvalue of tearing modes from MHD theory 
and found that the real part is consistent with the growth rate for tearing modes in the classical theory 
(Furth et al ., 1963). The imaginary part (or the frequency) is just in the range of Alfven waves (the 
phase velocity of tearing modes has die same magnitude as Alfven velocity). Therefore, Huang (1988) 



considered the plasma resistivity and the electric current along magnetic field line and derived the 
growth rate of Alfven waves, which is directly proportional to the plasma resistivity and the electric 
current density. The modulation of Alfven waves for the motion of a single electron possibly results 
in the resonance wave-partide interaction (Huang, 1989), and the growth rate of Alfven waves in this 
case is one or two orders of magnitude higher than the linear exponential growth rate and leads to a 
turbulent spectrum, which is a possible explanation for the anomalous growth of tearing modes. 

The following radio observations designed for coronal activity studies of the Sun and of solar-type 
coronae are in progress : 

a. Solar Observations 

In conjunction with T. Bastian at NRAO we are in the process of calibrating a unique database 
on two solar active regions. The data consist of VLA and OVRO observations. The OVRO observa- 
tions - made over a very wide range of frequencies - are complementary to the high spatial resolution 
VLA observations that were made at 10 frequencies in a 500 MHz band at 5 GHz and at 1.4 and 8 
GHz. This database will allow us to study for the first time the structure of active regions as a func- 
tion of frequency, and allow us to build a consistent model of the radio emission. However, analysis 
has been delayed by the discovery during the calibration procedure that the VLA’s calibration of solar 
observations was in error. The error has been traced and corrected. 

In conjunction with J. Neff at GSFC, we have obtained a very deep SWP high dispersion spectra 
of the Lyman-a line from a G type subgiant star. The unique aspect of this work is that the star has a 
radial velocity sufficiently high to remove the Lyman-a line core from the ISM absorption. Thus, we 
will be able to model the line emission from a star of spectral type similar to the Sun, but with a dif- 
ferent effective gravity. This procedure should allow us to understand more clearly the formation 
mechanism of this line in the Solar atmosphere. 


b. Stellar Observations 
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In collaboration with observers at Haystack, we initiated and performed a very sensitive VLB I 
experiment in an attempt to directly measure the sizes of the radio-emitting regions of stellar coronae. 
Of the six stars observed, three were detected and two were resolved. For two stars, we were able to 
place a limit on the size of the quiescent radio emission less than 0.4 stellar diameters. These are the 
first measurements of the sizes of stellar coronae. The success of these observations has also allowed 
us to immediately eliminate at least one possible radio emission mechanism (gyroresonance emission) 
from consideration. 

In conjunction with T. B^stian at NRAO, we are in the process of finalizing the analysis of 
several exceptionally large stellar radio flares on the M Dwarf AD Leo. These flares were remarkable 
in that the emission - despite the fact that it MUST have been produced by a coherent mechanism 
such as an electron-cyclotron maser or plasma radiation — was broadband in nature (8v/v > 0.03). 
This observation has forced a revision in the idea that the coherent mechanisms must be intrinsically 
narrow-band. As a consequence, it appears that the quiescent emission from M dwarfs may be attri- 
butable to a superposition of many small flare-like events. 
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disk He I X10S30 spectrohellograms have been available since Skylab from the National Solar Obser- 
vatory (NSO) at Kitt Peak. The Helium daU have been used, among other things, for studying the 
small regions which are called "dark points" when seen in Helium, under the working hypothesis that 
there Is a one-to-one association between dark points and XBP. This assumption Is based on a study 
carried out in 1973 (Harvey et at. 1975) using some of the earliest He I data available. This study was 
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Comments on the Observability of Coronal Variations 


L. Golub*, T. W. Hartquist^ and A. C. Quillen 1 ’^ 

1. Smithsonian Astrophysical Observatory, Cambridge Mass. 

2. Max-Plank-Institut fur Physik und Astrophysik, D-8046 Garching 

3. Astronomy Dept., Caltech, Pasadena Cal. 


ABSTRACT 

We discuss the observable variability of spectral lines in the soft x-ray and XUV region. 
Rapid variability of coronal emission, both in flaring and non-flaring structures has been reported 
and is particularly prominent when high spatial resolution is available. Examination of the ioniza- 
tion and recombination timescales for the formation and removal of ions with prominent Solar 
emission lines shows that variability timescales are often limited by these atomic processes, 
independent of the physical process which is causing the change in the Solar atmosphere. Rapid 
heating can lead to an initial freezing-in of abundances of some ions; observations of at least one 
low and one high excitation line from such an ion would permit studies of the time evolution of 
the emission measure and temperature. In some cases, rapid cooling leads to freezing-in of the 
abundance of an ion and observations of low excitation line of this ion will not yield accurate 
information about the thermal evolution. Thus, future observations of Mg X 609 A should be aug- 
mented by simultaneous observation at another wavelength, such as 63 A. In addition, with the 
ability to produce images in isolated spectral lines it becomes possible to select those for which 
rapid variability is observable, such as O VII, rather than lines which were selected on the basis 
of previous hardware constraints, such as O VIII. 


Submitted to Solar Physics 
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1. Introduction 

Techniques for the fabrication of thin film multilayer coatings capable of reflecting x-rays have 
progressed rapidly in recent years (for a review, see Baibee 1985). It has now become possible to 
design and construct astronomical instruments which employ multilayer coated optical elements and 
operate at normal incidence at wavelengths in the region from 44 to -304 A (Golub, et al. 1985, 
Underwood et al. 1987; Walker et al. 1988). With the development of normal incidence x-ray and 
XUV telescopes, high spatial (-0.1 arcsec) resolution observations of individual spectral lines emitted 
by the Solar coronal and transition region plasmas are now possible. Designs which include dispersive 
elements are also being developed, as are designs using multilayers at grazing angles permitting their 
use at shorter wavelengths (see the review by Catura and Golub 1988). 

Because this new technology is available and is beginning to be used in Solar observations, the 
present time is opportune for addressing questions related to the observability of temporal variations in 
individual emission lines caused by the dynamical and physical evolution of the emitting regions. In 
the following, we will identify previously detected Solar emission lines which are suitable for the study 
of variability in plasmas with the specific combinations of pressure, temperature and spatial scales 
which arc found in the Solar atmosphere, and then we will explore the types of variability which can 
be observed in each of these lines. 

Rapid non-flare variability of coronal emission from individual loop structures seen in Fe XV at 
284 A was reported by Sheeley and Golub (1978). These data had 2 arcsec resolution, which permitted 
the individual loops in active regions to be resolved, whereas simultaneous observations in broadband 
x-rays, which had somewhat lower resolution, did not resolve these structures. Short time scale varia- 
bility was not observed in x-rays and this lack of variability was interpreted as due to an averaging 
over the emission from several loops in the active region (Nolte, Solodyna and Gerassimenko 1979), 
which smooths out the variability of the individual structures. 

These observations show that rapid variability does occur within active regions, and that spatial 
resolution of better than 2 arcseconds is needed for studies of these structures. However, our present 
study shows that the interpretation given by Nolte, et al. for the lack of variability in x-rays may not 
be correct. The x-ray emission which they observed was dominated by Fe XVII in the thin Be filter, or 
by a combination of Fe XVII, O VII and 0 VIII in the Polypropylene filter, with the latter possibly 
dominating. As we discuss in §3.3, the time scale for variability of the emission in 0 VIII at the tem- 
perature and density appropriate to typical coronal regions is several minutes. This is equal to or 
greater than the time scales observed in Fe XV, for which ionization equilibrium is established in a 
few tens of seconds. 

Thus, the failure to observe variability in x-rays may have been an artifact of the instrumentation 
used, rather than anything intrinsic to the active region under study. Instruments which can provide 
imaging in a single spectral line are in some ways more susceptible to this problem, because an 
incorrect choice cannot be compensated by another line in the passband. For the study described above, 
the ability to choose O VII rather than O VIII would have provided diagnostics with 0.1 second time 
resolution. 

Observations in more than one spectral line are needed in order to distinguish between tempera- 
ture and density variations, since the emission from any single line by itself will yield information that 
mixes emission measure and temperature. For example, it is possible to observe temporal variability 
more rapid than the ionization times if the emission measure changes rapidly. 

High spatial and moderate spectral resolution observations may reveal not only the variability of 
individual loop structures, but also the possible substructure within coronal loops. In coronal heating 
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models based on reconnection and current dissipation (e.g. Kupems, Ionson and Spicer 1981), the pri- 
mary heating is confined to very thin current sheaths which may heat the loops in either a steady or 
intermittent fashion. In either case, the loops should exhibit small scale thermal substructure with 
associated variability. 

Martens et al. (1985) have reported observations of the corona using the HXIS instrument on the 
SMM satellite; their interpretation of the data is that "0.1% of the loop is filled with a plasma with a 
temperature of about 10^ K and a density of 10 1 * cm*^". The temperature and density in the remain- 
ing 99.9% of the loop volume are "much smaller". Observations such as these offer indirect evidence 
that loop substructure is present, and it is likely to be highly variable on short timescales. 

The aim of this paper is to discuss the appropriateness of specific XUV and soft x-ray emission 
lines for studies of coronal variability. Among our major conclusions are the realization that, for many 
of the typical regions found in the Solar outer atmosphere, the observable timescales for variability arc 
determined by the ionization and recombination timescales of the plasma rather than by the inherent 
nature of the physical process causing the variability. Moreover, we find that cases exist in which con- 
ductive cooling dominates and that in some of these cases the "freezing-in" of the ionization structure 
of a species affects the emissivity. If the particular lines being observed originates from a lowly excited 
level then the freezing-in will lead to the line’s strength varying slowly compared to the cooling time. 
On the other hand, we show that diagnostics exist which can be used to determine whether or not con- 
ductive cooling is important in a particular event. 

The organization of this paper is as follows: In section II we list the lines which are most likely 
to be used in future observations and describe their applicability in diagnostic studies of coronal emis- 
sion under the variety of conditions encountered on the Sun. In section III we consider the ways in 
which the timescale for electron impact production of particular species affects the utility of those 
species for diagnosing changes in the coronal plasma; in the same section, we also treat cooling plas- 
mas. 



-4- 


2. Selection of Lines 

The x-ray and XUV emission lines which we have studied were selected on the basis of several 
criteria: 

i) We chose the strongest lines in the Solar spectrum, based both on theoretical calculations and 
on observations which ensure that the line is actually present; 

ii) We arbitrarily chose to limit our search to the range of wavelengths for which multilayer 
optics might be used. We take this range to be ~15 A at the short end and 304 A at the other end. 

iii) Specific consideration was given to selecting those lines which would provide unambiguous 
diagnostics over the full range of temperatures from the chromosphere through flare levels (i.e. 10 to 
10 7 - 5 K). 

The results of this search are not entirely satisfactory, in that the combination of technological 
difficulties and placement in the spectrum of the most interesting lines overlap in unfortunate ways. 
The results of the search can be summarized by the following three groupings: 

I. 15*44 A: This region contains the emission lines which for pure coronal imaging and diagnos- 
tics might be considered the most desireable, viz. Fe XVII, 0 VII, 0 VIII, C VI, etc. However, it 
docs not seem possible at the present time to manufacture multilayers which work with any reasonable 
efficiency at normal incidence at these wavelengths. Grazing incidence (- 10 - 20°) applications arc 
possible, but these to some extent remove the enormous advantage in terms of spatial resolution which 
is achieved by going to normal incidence. 

II. 44-170 A: In this region production of high performance multilayers becomes feasible, espe- 
cially longward of 60 A. However, the coronal emission lines in this region arc comparatively weak 
and construction of a practical observing instrument is difficult because of the low intensity levels. 

III. 171 • 304 A: Strong emission lines with important diagnostic value are found in this spectral 
region and high reflectivity multilayers are easily manufactured at these wavelengths. However, the 
spectrum is relatively crowded with strong emission lines and, because only a few layers arc needed in 
order to achieve high reflectivity in the multilayer it is quite difficult to obtain a narrow band coating. 
Thus, the data will in general be contaminated by a mixture of different lines all imaged on top of each 
other. In principle this problem can be overcome by introducing a dispersive element into the optical 
system, or by very clever design of the multilayer, neither of these two solutions has yet been demon- 
strated in practice, however. 

The lines which we have chosen are listed in Table I-a. Sources examined in the selection pro- 
cess include, among others: Freeman and Jones (1970), Manson (1972), Malinovsky and Heroux 
(1973), Doschek and Cowan (1984) and Acton, et al. (1985). The decision as to which lines will be 
most useful for coronal studies depends to a large extent upon the plasma conditions which arc pre- 
valent in the atmospheric feature (or event) of interest. There is no single best answer to this question, 
since typical features or events usually have a connection to the photosphere or below, thus encom- 
passing temperatures from -5000 K to 10^ or as high as 10** K in extreme cases. However, a useful 
method of characterizing the study is to examine the coronal temperature, pressure and geometry, after 
which the lower temperature portions of the atmosphere are to some extent specified. 

A simplified tabulation of typical plasma parameters is provided in Table I-b. It should be under- 
stood that these are only typical values, listed in order to show the general range observed. In the 
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following discussion we use this table as a starting point for an examination of the spectral lines which 
will be most useful for studying each of these coronal phenomena. 

The Solar atmosphere is inherently complicated and atmospheric studies typically will require 
high resolution two-dimensional imaging and some type of spectroscopy. Structures with sizes ranging 
from the resolution limit (- 200 km) to several Solar radii (>10^ km) arc seen and major changes in 
these structures occur on times scales from milliseconds during the rise phase of a flare (bursts), to 
tens of minutes for flares (gradual phase), to months for the evolution of the large scale structure. It is 
not likely that a single instrument can be designed which will diagnose the full range of conditions 
which are present on the Sun, and in fact the usual answer is to provide a complement of instruments 
which together form a Solar observatory. 

The range of wavelengths over which multilayer coated optics can be used helps in deciding 
which subset of the overall problem will be attacked. In the "multilayer regime" the emission lines 
listed in table I-a cover a temperature range which includes the transition region and the corona, as 
well as thermal flare plasma. Thus our discussion centers around the type of problem typically treated 
by a grazing incidence x-ray telescope, with the addition of some of the lower temperature regions usu- 
ally treated by an XUV or UV spectroheliograph, such as the S-082A instrument on Skylab. The spec- 
tral lines which are of interest for the regions or events listed in Table I-b are obtained from Table I-a 
by comparing the temperature, emission measure and variability timescales; results are listed in Table 
1-b in the last column. 
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3. Response of the Plasma to Rapid Heating or Cooling 
3.1 General Principles 

Shapiro and Moore (1977) and Shapiro and Knight (1978) have pointed out that the plasma pro- 
perties in solar flare regions vary on timescales which are short relative to the timescales on which the 
ionization structure attains equilibrium. It is feasible to perform detailed calculations of the nonequili- 
brium ionization structure and emissivity for a small number of specific combinations of plasma condi- 
tions and thermal histories as did Shapiro and Moore (1977). 

However, coronal events display a wide variety of properties and in practice more general con- 
siderations, of utility for the study of any given coronal plasma, are required. For instance, Shapiro 
and Knight (1978) discussed the initial increases in the strengths of high excitation lines (i.e., those for 
which the excitation energy is comparable to or higher than kT^, where T Q is the original temperature), 
and the subsequent decreases due to collisionally induced ionization. They argued that the observed 
enhancement of a single properly-chosen high excitation line can be used to determine the increase in 
temperature, while the subsequent decay of the line strength yields information about the density of the 
plasma. They envisaged the ionization remaining frozen-in initially with its preheating structure and 
the increase in a line’s emissivity being due to the temperature increase; the contribution to the line 
strength’s increase due to stripping of less ionized species was ignored. The subsequent decrease in the 
line’s emissivity occurs at the density-dependent ionization rate. 

In fact, a plasma in which both the density and temperature have grown rapidly can be studied 
best by comparing the temporal evolution of the strength of a high excitation line and a low excitation 
line. The increase in the strength of the low excitation line will be due primarily to the increase in the 
emission measure (n g v, where V is the volume) since the temperature dependence of an electron 
impact excitation rate generally is fairly weak when the energy of the impacting electron is much 
greater than the excitation energy. The increase of the strength of the high excitation line will result 
from both the rise in the emission measure and from the rise of the temperature, the latter usually lead- 
ing to a rapid growth in the collisional excitation rate of a high excitation line. 

In structures in which the conductive cooling is very strong, the cooling may also occur on a 
timescale which is short compared to that required for electron impact induced ionization to produce 
species which normally would be abundant in an equilibrium situation. The absence of emission 
features which would be prominent in an equilibrium plasma at its temperature maximum may lead one 
to underestimate the maximum temperature attained. Qeariy, the variability of the plasma properties 
cannot be observed in lines of those species which cannot be ionized on timescales shorter than the 
evolutionary timescales of the plasma. 

Simultaneous observations of a high excitation line and of a low excitation line are also useful for 
diagnostic studies of a plasma which is cooling rapidly. In plasma which is cooling due to rapid con- 
ductive heat loss, the ionization states may be frozen in with a structure similar to that which obtains 
in a high temperature equilibrium plasma. If freezing in does occur, low excitation lines emitted from 
highly ionized species will continue to be strong and would generally increase in strength as a function 
of time if the plasma pressure remains constant during the cooling phase. However, the strengths of 
the corresponding high excitation lines would decrease relatively rapidly with plasma temperature. 
From simultaneous measurements of a low excitation line and a high excitation line of frozen in highly 
ionized species, the temperature and emission measure evolutions could be derived. As a practical 
matter, it is preferable if the two lines are of the same ion, in order to keep the analysis as uncompli- 
cated as possible. 
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An example of this effect is given in Figure 1, which compares the 52.9 A emission of Fe XV to 
the 284 A emission, using values taken from Mewe et al. (1985). The emissivity curves for the ion 
are identical, but dividing by the fractional abundance / for the respective lines in question introduces a 
major change, as shown by the upper curves in the figure. Thus, if material which is initially at a tem- 
perature of log T = 6.3 cools rapidly to a temperature of, e.g., log T = 6.0, then the decrease in 
strength of Fe XV at 284 A can be as much as a factor of three less than that of the 52.9 A line. 
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3.2 Variability 't imescales for Line Emission 


In section 2. we identified lines which have been detected previously and which are prime candi- 
dates for use in higher resolution studies of the corona. So that we can identify which lines are likely 
to be of diagnostic value, we have calculated the ionization and recombination timescales of various 
ions at different temperatures for plasma in which the pressure n £ T is arbitrarily taken to be equal to 
the value 3x10^ cm"^ K, where n g and T are the electron density and temperature. Data for the cal- 
culation of these timescales come from the compilation by Amaud and Rothenflug (1985) and refer- 
ences therein, primarily the papers of Woods, Shull and Sarazin (1981) and of Shull and van Steenburg 
(1982). In general, the time scales are inversely proportional to the pressure, but at high densities, 
diclcctronic recombination may be suppressed. 

The fraction, /, of an element that is in a given ionization stage in a low density equilibrium 
plasma can be found in the paper by Amaud and Rothenflug (1985). For several species, we give/ and 
various ionization and recombination timescales as a function of temperature in Tablc^2. In addition to 
several ions which are especially prominent in the Solar spectrum, the species C + and Fc +y are 
included because they have lines which arise from states with very different excitation energies, mak> 
ing them useful as diagnostics when studying rapid heating events. The species Fe +1 , O +0 and Mg +y 
arc abundant in equilibrium plasmas over a wide range of temperatures and are observable in at least 
two types of Solar coronal regions (viz. Table 5). 


We calculate timescales for four processes which can influence the appearance and disappearance 
of a line. The timescale, tjp for the ionization of the next less ionized stage of an element and the 
timescale, tj r , for the ionization to the next more ionized stage of that element are shown, t^ is the 
timescale for the recombination to the next less highly ionized stage of the element and t^ is the 
timescale for the recombination of the next more highly ionized stage of the element. This nomencla- 
ture is defined schematically in Figure 2. The timescales are normalized at n fi T = 3x10 5 cm’ -K. 
The ionization timescale is equal to (n e [C ea (T) + C^CT)])' 1 where C ga is the autoionization excitation 
rate coefficiant (significantly large for the lithium like ion 0+5 and for iron ions in the Sodium - Sul- 
phur series) and is the direct (electron impact) ionization rate coefficient The recombination 
timescale is equal to (n e (o r (T) + a^CT)])* 1 where a r is the radiative recombination rate and is the 
dielectronic recombination rate. Table 3 contains results extracted from Tables 2 for the temperature 
T max at which a speeds attains its maximum abundance. 

Normally the power emitted per unit volume in a single line by a plasma in equilibrium is written 
as the product of an emission rate coefficient, F\ , and n g . Note that even though in equilibrium plas- 
mas the concentrations of species such as C*® and Fe +9 may peak at the same temperature, the 
behavior of the line emissivities will differ greatly in response to rapid heating, as shown by the differ- 
ences in the P yjf curves; this behavior during rapid heating was discussed at the beginning of §3.1. 
?ylf should be multiplied by the fraction of the element in the appropriate species to calculate the 
power in the line when the plasma is not in equilibrium. 


Assume that we are interested in observing the emission from a species formed by the ionization 
of more lowly ionized species in a plasma which has been rapidly heated. The charge on the nucleus 
of that species is z and its ionization stage is £ where £ = 0 for a neutral species and £ = z for a fully 
ionized system. The number density is written n z ^. Because the ionization timescales increases 
rapidly with charge, it is normally sufficient to assume that a large fraction of an element is initially in 
the next lowest ionization stage when calculating the evolution of the population of a species. During 
the interval when n r peaks, the system of equations 

dn z.c-l /dt = “ n z,C-l A if + n z.C A rr 
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dn z,C /dt = n z,C-l A if + n z.C + l A rf 
-< 1A *ir +1/ 0 n z,C 


n z,C+l 


• n z" n 


'z,C " V-i 


n z,C-l (t=0) = n z 

^ jft 4c 

will govern the behavior of n z The asterisks on t .p t ^ , t - r . and t n indicate that they should 
be evaluated by multiplying the results in Table 2 by ( 3x10 cm’ -K / n g T ). Values of P^/f can be 
used with the derived n r to calculate the line strengths. 

The abundance of a species n z r in a plasma which is heated can be approximated from the 
above equations if a few simplifying assumptions are made. We assume that initially n z ^ = n z , and 
that throughout the initial rise in n z * that n z j can be set to 0. The question we ask is, ifthe plasma is 
heated instantaneously to a higher temperature, what is the timescale to form species z,£? Given these 
assumptions, one can show that 

"z,C = (n z T app ' 1 if* ( 1 " ® ^ 

with 


Spp = 


+ t, 


-1 


rr 


Values of this quantity for the species which we considered are given in Table 3. 

The abundance of a particular species in a recombining plasma cannot always be followed with 
such a small set of equations, because the recombination time scale is not a strongly monotonically 
decreasing function of £. At log T = 6.5 the values of for Fe + ^"\ Fe + ^, Fe + , and Fe + ^ arc 
43.4s, 27.8s, 14.5s, and 21.8s respectively. Hence, in a plasma which has cooled very rapidly, the 
recombination eventually leads to significant populations in a range of ionization stages. However, in 
many cases it is sufficient to compare the recombination timescales to the cooling timescales to deter- 
mine at what temperature the ionization structure is last in statistical equilibrium. As long as the 
recombination continues to occur on a timescale long compared to either the cooling timescale or the 
timescale since rapid cooling has stopped, the ionization structure can be assumed to be that which 
obtains in an equilibrium plasma at the temperature at which cooling first becomes more rapid than 
recombination. 

The set of equations given above should be sufficient to follow the abundance of a species z,£ in 
a cooling plasma under cenain restricted conditions. For instance, if a plasma is cooling from an equili- 
brium at which n z * » n z for all i>2, the three equations should give a reasonable description of 
the evolution of n z t. Clearly, the initial values of n z r_j> n z r and n z j should be specified to be those 
for a plasma in equilibrium at the initial temperature. ’ 

We have given in Table 4 the timescales for cooling due to radiative losses in equilibrium plasma 
in which n g T = 3xl0 15 cm -K initially; these timescales are appropriate for a plasma cooling isochor- 
ically and were obtained by using data given in a figure produced by R. Edgar with J. Raymond’s code 
and provided by J. Raymond, showing the power loss coefficient as a function of temperature. The 
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cooling timescale for isobaric cooling is obtained by multiplying the results given in Table 4 by 5/3. 


The timescale for conductive cooling may be calculated by considering a loop of half-length 1 in 
which radiative losses are balanced exactly by heating and which is cooling isochorically due to con- 
duction. Several approximate solutions have been calculated for this case; Krall has suggested that the 
temperature, T , at the apex of the loop evolves as 


T = T 

A o A oo 


(1 


1.4X10 9 T* 2 
nl 2 ' 


where is the initial temperature at the apex of the loop. 

In calculating the conductive cooling time we will use Krall’s result. However, we note that the 
equation of thermal balance is a partial differential equation with time and spatial derivatives; if one 
assumes that conductivity is the only process responsible for energy loss or gain, that flows are very 
subsonic, that the loop is a cylinder in which conductive heat transfer is parallel to the axis of sym- 
metry, and that the gas cools either isobarically or isochorically, then the equations can be solved by a 
separation of variables. No solution exists for a cooling loop with T = 0 at the base and T a maximum 
at the top. Consequently, any value for the conductive cooling time which one may adopt is somewhat 
arbitrary and is likely to be corret only to within an order of magnitude of the true value. 
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3.3 Emission From Solar Coronal Features 

Wc have given in Table 5 typical timescales for the appearance and disappearance of lines prom- 
inent in the spectra of the coronal features listed in Table lb, using the temperature and density values 
appropriate to the different types of regions under consideration. The timescales given take into 
account the number, n z ^ of atoms contributing to the emission and arc calculated in several steps, 
x , the minimum time for the appearance of the line, is the ionization time for formation of -60% 
(one e-folding) of the ion from the next less ionized state of the element, calculated assuming the tem- 
perature instantaneously becomes the maximum temperature for emission of the line. These times 
were calculated from the times listed in Table 3 using the temperatures and pressures typically found in 
the coronal features emitting the line, x recom is a lower limit on the disappearance time of the line. 
It is the time for 60% of the ion to recombine to the next less ionized state at the temperature and 
pressure found in the coronal feature. 

For comparison also given is x ^ CQn( j , the disappearance timescale of the line assuming that 
conductive cooling is the dominant cooling mechanism, and x ^ ra( j . the disappearance timescale of 
the line assuming that radiative cooling is the dominant cooling mechanism. These are the times for 
the emissivity of the line to drop to 1/e of its maximum value. Conductive cooling rates were calcu- 
lated assuming constant heat flux throughout the loops (Krall 1977). Radiative cooling rates were cal- 
culated using analytic fits to Raymond’s power loss function (viz. Rosner, Tucker and Vaiana 1978). 

T dis cond 311(1 T dis rad were ca ^ cu ^ ate< ^ assuming ionization equilibrium. In these cases, it is 
necessary to consider not only the timescale for a given change in temperature which is given directly 
from the cooling curves, but also to consider the change in emission which occcurs for a given change 
in temperature. In this way, the timescales given in Table 5 for disappearance of a line refer to the 
same quantities listed for appearance of the line. 

The cooling timescale due to radiative losses is inversely proportional to the pressure, in optically 
thin plasmas. For the high temperature plasmas encountered in active regions and flares in the Solar 
corona, the pressures are usually very much higher than the value of 3x10 13 which we have used as a 
normalization. Thus radiative cooling times for flares are typically 10^ to 10^ seconds and arc compar- 
able to the typical conductive cooling times calculated for flare loops. We wish to point out only that 
conductive cooling may in some cases be more important than radiative cooling in the coronal plasma. 
As discussed above, freezing-in can occur when the cooling is very rapid, and in those cases simul- 
taneous observations of properly chosen high and low excitation lines should be used in order to diag- 
nose both the emission measure and the temperature variations. 
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3.4 The Case of Mg X 

Of the lines listed in Table 5, the only one for which freezing-in might occur under coronal con- 
ditions is Mg X emitted by large scale coronal structures. However, as pointed out in §3.3, the correct 
values for conductive cooling times are uncertain by a large factor and conductive cooling might be 
important in a wider variety of situations in the Solar atmosphere; the opposite may also be true. Under 
the assumption that there will turn out to be cases for which we need to consider rapid cooling we give 
a treatment which is of general validity for a gas which has just begun to cool isochorically. 

The total power emitted in a line in a cooling plasma is given by 

Ex(z,Q = Px(z,Q X(z,0 (n. 2 V) 

where P^(z,Q is the line power rate coefficient, V is the volume of the emitting gas, X is the 
wavelength, and 

X(z,0 = n^/n. 


In general, the mass of gas will remain constant as it cools, i.e., n e V remains constant. We restrict con- 
sideration to rapid cooling so that n g remains constant; this is equivalent to requiring that the sound 
crossing time in the loop is long compared to the cooling time. The timescale on which the line 
strength varies is then given by 

l dis = ( P *‘ Iff dT/dt + X ~' dX/dt >"' 

A reasonable estimate for t djs is obtained by replacing dT/dt with 

( T/t^ + T/Wa* ) 


and X’ 1 


dX/dt with the minimum of 


[ 


1 dX^z.Q T 
Xeq(z,Q dT tnd^ool 


tcond,cool 



where X is given from equilibrium calculations. 

The data necessary to calculate t d j s for Mg X are given in Table 2 and Table 4, and by the equa- 
tion for the conductive cooling time. If the Mg + ^ abundance freezes-in then observations of Mg X at, 
for example, both 63 A and 609 A will permit a determination of the temperature over a broader range 
than if the Mg attained its equilibrium value. This also means that the timescale for disappearance 
of the Mg X emission can be significantly larger than that of another ion formed at a similar tempera- 
ture. Under conditions of rapid cooling (which presumably alternate with periods of heating) the 609 
A line will show less variability than will the 63 A line. Measurement of both lines permits an accurate 
temperature determination, and measurement of a line of a species which attains equilibrium abundance 
rapidly (such as O VII) can be used to determine the emission measure and thence n . 

V 
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4. Summary 


With a view toward observations of the Solar transition region and corona using imaging of indi- 
vidual spectral lines in the soft x-ray and XUV portion of the spectrum, we have discussed some of the 
major factors which enter into the choice of spectral line to be used. The discussion has centered 
around the lines which are available to the technology of multilayer coatings, but is more general in its 
range of applicability. Because of our interest in coronal heating and the rapid variability of coronal 
structures we have concentrated on the observability of variations in line strength for the range of 
features found in the Solar corona. Examination of the ionization and recombination timescales for the 
formation and removal of ions which are prominent in the coronal spectrum shows that variability 
timescales are often limited by these atomic processes, independent of the physical process which 
might be assumed to be causing the changes. 

In some situations, the plasma may be heated on timescales which are short compared to the time 
needed for the ionization structure to attain equilibrium. Such rapid heating can lead to a freezing-in of 
the initial ionization structure. In such cases, comparison of a low excitation line of a species to a high 
excitation line (preferably of the same species) would permit accurate determination of the time evolu- 
tion of both the emission measure and the temperature. Similarly, there are some situations in which 
rapid cooling may lead to a freezing-in of the abundance of an ion. In particular. Mg X viewed in 
large scale coronal structures should show this effect. Observations in the 609 A line of Mg X should 
be supplemented by the short wavelength line such as Mg X at 63 A. 


Studies of rapid variability in active regions should be performed' with spectral lines which are 
capable of being formed on short timescales. For example, at typical active region temperatures of 2- 
3X10^ K, the species O VII can appear on times a factor 100 shorter than those for the appearance of 
0 VIII. Disappearance times for 0 VII are about one order of magnitude shorter than those for O VIII. 


This work was supported by NASA Grant NAGW-112 to the Smithsonian Institution and by the 
Smithsonian Institution Visitors’ Program. We would like to thank Dr. John Raymond for helpful dis- 
cussions during preparation of this manuscript. 
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6. Figure Captions 
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Figure 1. Comparison of the emissivity of a low and high excitation line of the species Fe 
The upper curves show the emissivity divided by the fractional abundance for each line. 


Figure 2. Schematic illustration of the nomenclature used for labelling the formation and 
removal processes of species i. 





Table I-a. Soft X-ray and XUV Lines Included in This Study 


Species Wavelength (A) l°g(T max ) 


Fe XVII 

16.78 

6.6 


17.05 



17.10 


0 VIII 

18.97 

6.5 

0 VII 

21.60 

6.3 


21.81 



22.10 


C VI 

33.74 

6.1 

Si XII 

45.52 

6.3 


45.68 


Si XI 

49.22 

6.2 

Mg X 

57.9 

6.1 

Si VIII 

61.0 

6.0 

Fc XVI 

62.88 

6.4 


63.71 


Mg X 

63.15 

6.1 


63.29 


Nc VIII 

88.08 

5.8 

Fc XVIII 

93.93 

6.8 

Nc VIII 

98.26 

5.8 

Ne VII 

115.52 

5.7 


116.69 


Fe XXII 

116.29 

7.1 


117.17 


Fc IX 

171.07 

6.0 

0 V 

172.17 

5.4 

0 VI 

172.94 

5.5 


173.08 


Fc X 

174.53 

6.1 


177.24 


Fe XI 

180.41 

6.2 


182.17 



188.22 


0 VI 

183.95 

5.5 


184.13 


Fe XXIV 

192.03 

7.3 

Fe XII 

192.40 

6.2 


193.52 



195.13 


Fe XIV 

211.32 

6.2 


274.24 


He II 

256.37 

4.6 

Fe XV 

284.16 

6.3 

He II 

303.79 

4.6 

Si XI 

303.33 

6.2 


Tabic I-b. Typical Plasma Parameters for Solar Atmospheric Features 


Target 

Size [xl0 9 cm] 

cn 

i 

i 

T [10 6 K] 

Lines 

Quiet Corona 

20 

l-3xl0 8 

1.6 

0 VII, C VI, Mg X. 
Fe XI, Fe XII 

Extended Corona 

70 

10 7 

1.4 

Fe X, Mg X, Ne VII 

Polar Plumes 

50 

10 8 

1.3 

Mg X, Fe X, Ne VIII 

Active Region 

5-20 

10 9 -10 10 

2-3 

0 VIII, Fe XVII, 
Fe XVI. Fe XV 

Loop Microstructurc 

<1 

10 11 

10 

Fe XVIII. Fe XXII, 
Fe XXIV 

X-ray Bright Point 

.2-2 

5xl0 9 

1.8 

0 VII, Si XII, 
Fe XIV 

Compact Flare 

2-3 

10" 

20 

Fe XXII, Fe XXIV 

Two-Ribbon Flare 

10 

2xl0 10 

10 

Fe XVIII, Fe XXII, 
He II 



Table 2. Ionization and Recombination Timescales as a Function of Temperature 


log T 

l if 

l ir 

l rr 

l rf 

-log/ 



C VI 




5.7 

7.8e3 

2.1e5 

1.6e2 

86 

1.7 

5.8 

1.4e3 

2.4e4 

2.2e2 

1.3e2 

0.9 

5.9 

3.5e2 

4.2e3 

2.6e2 

1.9e2 

0.4 

6.0 

1.2e2 

l.lc3 

2.7e2 

2.9e2 

0.2 

6.1 

56 

4.0e2 

2.7e2 

4.3e2 

0.4 

6.2 

31 

1.8e2 

2.8e2 

6.5e2 

0.7 

6.3 

20 

1.0e2 

3.1e2 

9.9e2 

1.0 

6.4 

15 

65 

3.7e2 

1.5e3 

1.4 

6.5 

12 

48 

4. 7e2 

2.3e3 

1.7 



Fe IX 




5.6 

95 

69 

0.76 

0.37 

2.2 

5.7 

4.0 

19 

0.89 

0.51 

1.3 

5.8 

2.0 

7.2 

1.1 

0.70 

0.8 

5.9 

1.2 

3.4 

1.5 

0.96 

0.5 

6.0 

0.79 

2.0 

2.2 

1.4 

0.5 

6.1 

0.60 

1.3 

3.3 

2.0 

0.9 

6.2 

0.50 

0.99 

5.1 

2.9 

1.8 

6.3 

0.45 

0.82 

8.1 

4.5 

3.5 



Fe XVIII 



6.3 

5.5e3 

1.2e4 

l.le2 

41 

2.5 

6.4 

1.4e3 

2.8e3 

1.3c2 

59 

1.5 

6.5 

4.8e2 

8.9e2 

1.5e2 

82 

0.9 

6.6 

2.2e2 

3.7e2 

1.7e2 

l.le2 

0.5 

6.7 

1.2e2 

1.9e2 

2.1e2 

1.5e2 

0.4 

6.8 

77 

1.2e2 

2.8e2 

2.0c2 

0.6 

6.9 

57 

85 

3.9e2 

2.8e2 

1.0 

7.0 

47 

68 

5.7e2 

4.1e2 

1.8 

7.1 

42 

59 

8.6e2 

6.1e2 

3.0 



Mg X 




5.8 

2.7e2 

1.3e3 

13 

58 

3.0 

5.9 

90 

3.7e2 

21 

88 

1.5 

6.0 

39 

1.4e2 

34 

1.3e2 

0.7 

6.1 

21 

70 

56 

2.0e2 

0.7 

6.2 

13 

41 

92 

3.1e2 

0.9 

6.3 

9.6 

28 

1.5e2 

4.6e2 

1.2 

6.4 

7.8 

22 

2.5e2 

6.3e2 

1.5 

6.5 

6.8 

18 

4.1e2 

7.8e2 

1.7 

6.6 

6.5 

17 

6.7e2 

9.0e2 

1.8 

6.7 

6.5 

16 

l.le3 

1.0e3 

1.9 

6.8 

6.9 

17 

1.8e3 

1.2e3 

2.2 




O VII 


5.3 

5.6e-2 

- 

5.4 

6.5C-2 

- 

5.5 

7.7e-2 

- 

5.6 

9.2e-2 

6.e9 

5.7 

0.11 

8.e7 

5.8 

0.14 

2.7e6 

5.9 

0.16 

1.8e5 

6.0 

0.20 

2.3e4 

6.1 

0.24 

4.5e3 

6.2 

0.29 

1.3e3 

6.3 

0.36 

4.8e2 

6.4 

0.44 

2.3e2 

6.5 

0.54 

1.3e2 

6.6 

0.68 

88 



0 VIII 

5.9 

1.8c5 

3.2e6 

6.0 

2.3c4 

2.6e5 

6.1 

4.5c3 

3.9c4 

6.2 

1.3e3 

8.6c3 

6.3 

4.8e2 

2.7e3 

6.4 

2.3c2 

l.le3 

6.5 

1.3c2 

5.5e2 

6.6 

88 

3.3e2 

6.7 

66 

2.3e2 

6.8 

54 

1.8c2 

6.9 

49 

1.5e2 

7.0 

47 

1.4c2 

7.1 

48 

1.4e2 

7.2 

50 

1.4e2 

7.3 

55 

1.5e2 


28 

11 

2.8 

42 

17 

1.1 

63 

26 

0.2 

94 

40 

0.1 

1.4e2 

61 

0.0 

2.1e2 

93 

0.0 

3.0e2 

1.4e2 

0.0 

4.0e2 

2.0c2 

0.0 

4.4e2 

2.7c2 

0.0 

4.3e2 

3.3c2 

0.1 

4.2e2 

3.9c2 

0.3 

4.4e2 

4.6c2 

0.6 

5.0e2 

5.7c2 

1.1 

6.2c2 

7.4c2 

1.7 


1.4c2 

94 

3.1 

2.0c2 

1.4c2 

2.1 

2.7c2 

2.1c2 

1.3 

3.3c2 

3.1c2 

0.7 

3.9c2 

4. 7e2 

0.4 

4.6c2 

7.1e2 

0.3 

5.7e2 

l.lc3 

0.5 

7.4c2 

1.6c3 

0.8 

1.0c3 

2.5c3 

1.1 

1.4c3 

3.8c3 

1.3 

2.1c3 

5.8c3 

1.6 

3.2c3 

8.9c3 

1.8 

5.0c3 

1.4e4 

2.0 

7.9c3 

2.1e4 

2.2 

1.3e4 

3.3e4 

2.3 




*) At the temperature of maximum emission. 






Table 4. Radiative Cooling Timescales for an Equilibrium Plasma 
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Tabic 5. Heating and Cooling Timescales for Solar Coronal Features. 


Feature 

Line 

\pp. 

T recom. 

T dis. cond. 

x dis. rad. 

Quiet Corona 

0 VII 

3.0 

3xl0 3 

5xl0 4 

10 4 

Extended Corona 

Fe X 

120 

500 

2xl0 4 

1.5x10 s 


Mg X 

2000 

4xl0 4 

2x!0 4 

1.5x10 S 

Polar Plumes 

Fe X 

12 

50 

10 5 

1.4xl0 4 


Mg X 

370 

4.5xl0 3 

10 5 

1.4xl0 4 

Active Region 

0 VIII 

28 

200 

8xl0 4 

I0 3 


Fe XVII 

4.0 

30 

10 5 

10 3 

Microstructurc 

Fe XVIII 

1.0 

1.2 

600 

600 


Fe XXII 

3.0 

0.7 

250 

400 

X-ray Bright Point 

0 VII 

0.1 

100 

800 

500 

Compact Flare 

Fe XXII 

3.0 

1.3 

400 

10 3 


Fe XXIV 

10 

2 

750 

2xl0 3 

Two-Ribbon Hare 

Fe XVIII 

5.0 

6 

10 4 

3xl0 3 


Fe XXII 

15 

3 

5xl0 3 

2xl0 3 


*) Using T and n values from Table I-b. All times are in seconds. 
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ABSTRACT 

We use spectra of the 7100 and 8860 A bands of TiO to estimate the temperature and projected 
area of cool spots on two active K dwarf stars: HD 82558 and the flare star BD +26° 730. 

OBSERVATIONS AND ANALYSIS 

The TiO bands are sensitive indicators of temperature in stellar atmospheres. The very presence of 
these features in dwarf stars with T e f/ > 4250 K («K5V) is direct evidence for cool spots. Ramsey and 
Nations (1980) flrst used the 8860 A TiO band to detect spots on HR 1099. Following their success, 
work has concentrated on RS CVn stars (Huenemoerder and Ramsey 1987; Huenemoerder 1988). In 
this paper we discuss the extension of this technique to measurements of the spot characteristics on 
active dwarf stars. 

We obtained spectra of the TiO region using the the National Solar Observatory’6 McMath Solar 
telescope and the stellar spectrograph system described by Smith and Jaksha (1984). Using the Bausch 
and Lomb grating in flrst order and the 105 mm transfer lens, the 2 pixel resolution on the TI CCD 
detector was A/AA=13,000. Signal- to-noise ratios were typically 100 to 300. 

The two wavelength regions we selected, 7100 and 8860 A, contain TiO absorption bands with 
various temperature sensitivities. The 7100 A bands (at 7055, 7088, and 7126 A; the 7 (0,0) system) 
flrst become visible for an immaculate 6tar at spectral type about K5, while the 8860 A band (the 
6 (0,0) system) is not evident until Ml. Modeling the two bands simultaneously, therefore, permits 
an estimate of both the size (from the absolute band strengths) and the temperature (from the ratio 
of the band strengths) of cool spots on the stellar surface. Our approach is essentially a simplified 
version of the one used by Huenemoerder and Ramsey (1987), in which the normalized spectrum from 
a spotted star ( F to t a i ) is modeled as the weighted sum of ^ tpot and Fphot, the normalized spectra of 
suitable standard stars with the effective temperature of the spot and the photosphere, respectively. 
The model is given by 


ri _ AR\F t tpot + (1 — A)Fphot 

Ftotal = arT+~i~A) ' 

where A is the total fractional projected area of spots on the observed hemisphere (weighted by limb- 
darkening) and R\ is the surface flux ratio between the spots and the photosphere. We assumed that 
R\ can be represented by the ratio of the respective Planck functions: i.e., R\ = Bx(T tpot )/B\(T p h,ot)- 
We performed simultaneous least squares fits in both wavelength regions to the above equation. Various 
standard star spectra were used to simulate F tpot and F p h ot until the best fit was achieved. Values of 
T e // for the standard stars were taken from de Jager and Nieuwenhuijzen (1987). 
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RESULTS AND DISCUSSION 


In this preliminary study, we have focussed on two dwarf stars: BD +26°730 (=Gliese 171.2A; 
K5Ve, B-V=1.12) and HD 82558 (K2Ve, B-V=0.91). We used the spectrum of f Boo B (K4V, B- 
V=1.16) as the model immaculate photosphere for BD +26°730 and the spectrum of HD 149661 (=12 
Oph; K1V, B-V=0.82) for the immaculate photosphere of HD 82558. A sequence of models were then 
generated for each star, using M giant spectra of varying T e / / to simulate spot spectra (the TiO bands 
are insensitive to gravity; Wing, Dean, and MacConnell 1976). 

Excess absorption in the TiO bands is clearly present in the spectrum of BD +26°730 (Fig. 1). The 
best fit to both bands was obtained using 2 Cen (M4.5ni) to simulate the cool spot spectrum. With 
T pho t = 4250 K, we find T. pot =3500 ± 150 K, AT s T phot -T tpo t=750 ± 150 K and A = 0.20 ± 0.04. 

The 7100 A band of HD 82558 is somewhat weaker than that of BD +26°730 (Fig. 2a), while the 
8860 A band does not appear to be present (Fig. 2b) at all, indicating that the spots on HD 82558 
must be warmer than those on BD +26°730. We found that fi Gem (M3III) was the coolest star whose 
spectrum we could use as F 4p<)t and still fit the HD 82558 data. With T p hot — 4900 K, we derive T, pot > 
3700 K, AT < 1200 K, and A >0.22. 

BD +26°730 (Bopp et al. 1983) and HD 82558 (Fekel et al. 1986) are among the most chromospher- 
ically active K dwarfs known. Both have rotational periods of less than 2 days. BD +26°730 exhibits 
flare activity, which is unusual for K dwarfs, and shows large, long-term luminosity changes, probably 
arising from a stellar activity cycle (Hartmann et al. 1981). The derived area coverages are similar to 
those found on the most active BY Dra and RS CVn variables (e.g., Poe and Eaton 1985; Rodono et 
al. 1986; Huenemoerder 1988). The derived temperature differences between 6pot and photosphere fall 
in the range seen for other BY Dra variables (e.g., AT = 600 K for BY Dra, Vogt 1983; AT = 850 K 
for AU Mic, Rodono et al. 1986). 

Spectroscopic determinations of spot properties measure absolute 6pot parameters, unlike photo- 
metric methods, which are differential in nature. We are continuing TiO observations, both in an effort 
to survey the spot properties of flare stars and to search for rotational modulation. Analysis of the 
latter will allow us to partly remove the geometric ambiguities in the interpretation of A and determine 
actual positions and sizes of the dominant starspots. With the survey results, we can explore (among 
other things) whether the increased area covered by starspots could explain the greater number and 
strength of stellar flares relative to the Sun. 

This work is based on data obtained at the National Solar Observatory (NSO), which is operated 
by AURA, Inc. under contract with the NSF. The research is supported by the Smithsonian Institution 
Fellowship Program, NASA grant NAGW-112, and a fellowship from the National Research Council. 
Finally, we thank the NSO for granting us observing time and travel support, and for their hospitality. 
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Figure 1. The 7100 A (a) and 8860 A (b) spectral regions of the K5Ve flare star BD +26°730 
(upper solid curves), the photospheric proxy { Boo B (K4V; dashed curves), and the starspot proxy 
2 Cen (M4.5III; lower solid curves). Note the excess absorption in the BD +26°730 spectra at the 
wavelengths of the TiO bands (which are clearly visible in the 2 Cen spectra). We derive T tpot = 3500 
± 150 K and A = 0.20 ± 0.04. 
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Figure 2. The 7100 A (a) and 8860 A (b) spectral regions of the K2Ve flare star HD 82558 (upper 
solid curves). The spectrum of the photospheric proxy HD 149661 (K1V) is shown in 2a (dashed line). 
Note the absence of TiO absorption at 8860 A, indicating warmer spots. An M3III star was the coolest 
that would fit the data at both wavelengths (/i Gem; lower curves), implying T apot > 3700 K, and 
A >0.22. 
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Abstract 

We study the current sheet created in the wake of an erupting filament during a two- 
ribbon flare. A comparison with the geomagnetic tail shows that the physics of these systems 
is very similar, and that the existence of super Dreicer fields and the generation of neutral 
beams travelling down the postflare loops with small pitch angles may be expected. The 
observational evidence for neutral beams is reviewed, and found to be generally supportive, 
while contradicting the widely held hypothesis of electron beams. 

Numerical simulation of individual particle orbits confirms the similarity with the ac- 
celeration process in the geomagnetic tail. It is further demonstrated that the results for 
selfconsistent numerical simulations of the current sheet in the geomagnetic tail can directly 
be scaled to the coronal current sheet, and the scaling parameters are derived. 

Subject Headings: Plasmas, Particle Acceleration, Sun: Magnetic Fields, Sun: Flares, Sun: 
X-Rays 

* Paper presented at the “Second Workshop on Impulsive Flares' 1 , University of New Hampshire. Sept. 19** 
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1 Introduction 


Two-ribbon flares are the largest flares observed on the sun, and filament eruptions are often 
considered as their physical cause (Carmichael 1964. Van Tend and Kuperus 1978). Observations 
(Dennis 1985, Oliver et al. 19S6), as well as numerical models (Kuin and Martens 1986, Mikic 
et al. 1988), and theoretical considerations (Pneumann 1980, Kaastra 1985) demonstrate the 
origin of a current sheet below the erupting filament, where the postflare loops are formed, and 
where most of the conversion of magnetic energy to particle kinetic energy takes place. This 
current sheet will be considered here in detail, and compared to the geomagnetic tail and a 
laboratory experiment (Stenzel et al. 19S5). 

During a large two- ribbon flare about 10 32 erg of magnetic energy is dissipated and eventually 
radiated away over the whole electromagnetic spectrum, or emitted as relativistic particles. The 
peak conversion rate, in the impulsive phase of the flare, is about 10 29 erg s -1 . In the current 
sheet the magnetic energy is either converted into thermal energy, producing a superhot (10 s K) 
plasma in the flare loops, or into electron beams, or. perhaps, neutral beams. The observed hard 
X-rays would be radiated directly from a superhot plasma, while electron beams would generate 
‘‘thick target" hard X-rays upon impact on the chromosphere, after traveling downward along 
the postflare loops. There has been quite a debate in the literature on the “thermal'’ (superhot 
plasma) versus the “nonthermaT (electron beam) hypothesis (see Dennis 1988), and neither 
theory is fully selfconsistent or fully consistent with observations at the moment. 

In this paper we will further explore the third alternative mentioned above: neutral beams. 
This was originally proposed by Simnett and Strong (1984) and Simnett (1986), because these 
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beams would explain several observations, such as the timing of the chromospheric upflows versus 
the hard X-rays, much better than either electron beams or a superhot plasma. Recently Martens 
(19S8a,b) and Martens and Young (1988) demonstrated an acceleration mechanism for these 
beams that is analogous to the acceleration mechanism that operates in the geomagnetic tail 
(Speiser 1965, 1967, 1968). Several drawbacks of pure electron beams, such as the unacceptably 
large currents and self magnetic fields induced by the beams (Hoyng et al. 1976), as well as the 
reflection or “mirroring” of the beam electrons with suffiently large pitch angles, are nonexistent 
in this model because the beams are neutral - with equal amounts of protons and electrons 
traveling at the same speed - and the beams are injected into the postflare loops with very 
small pitch angles, as they are in the geomagnetic tail (see Lyons and Williams 1984 for a recent 
review). 

In Section 2 of this paper we will present a brief summary of the acceleration mechanism 
for the neutral beams, as it is proposed for the geomagnetic tail and in the current sheet of 
two-ribbon flares. Essential for this mechanism is the existence of very strong electric fields 
(much larger than the Dreicer field) perpendicular to the magnetic field inside the sheet. W» 
will analyze the the theoretical and observational evidence for the existence of these fields. A 
comparison between the physical parameters of the coronal current sheet and of the geomagnetic 
tail will further reveal that despite the enormous differences in their physical parameters, the 
non-dimensional description is very similar, and consequently an identical acceleration proccess 
may indeed be expected. 

In Section 3 we will explore the observational consequences of neutral beams. Dennis (19SS) 

has recently compared these beams to the Loch Ness monster; because it is almost impossible 

,\ 

\ 
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to find observations that rule them out directly. We take issue here with this assertion by 
discussing three independent tests for neutral beams with protons in the 100 KeY range. The 
first test is the observation of redshifted Lymano radiation generated by transitions between 
protons within the beam and neutral chromospheric hydrogen, a test that was proposed by 
Orral and Zirker (1976) and Canfield and Chang (1985). This crucial test could be carried out 
on a UV spectrometer with modest temporal and spatial resolution (1%). The second test is 
the measurement of the polarization of line emission caused by the beam impact ionization of 
chromospheric hydrogen. This observation has been carried out recently at Meudon in the Ho 
line (Henoux and C'hambe 19S8, Henoux et al. 1988, and Henoux 1989) for three different flares. 
The results are consistent with the neutral beam model, and totally inconsistent u-ith de hi Ke V 
electron beams. The third one is the observation of impact polarization in the S I 1437 A line 
with the VYSP instrument aboard SMM (Henoux el al., 19S3), which leads to a similar result. 
In this Section we will also summarize the evidence for neutral beams offered by Stinnett 1 19>G'. 
and the measurement of an upper cutoff in proton spectra (Heristchi et al. 1976). expected in 
the present model. 

Finally, in Section 4, we will present our numerical results for the simulation of particle 
orbits in the current sheet of a two- ribbon flare. The results are consistent with the numerical 
Tesults of Speiser (1965, 1967) for particle orbits in the geomagnetic tail. We will then explore 
the requirements for fully selfconsistent models for the current sheet, that take into account 
the current and field generated by the accelerated particles. The familiar MHD-description is 
inadequate for this type of current sheet, since the particles are essentially collisionless during 
and after the acceleration, just as in the geomagnetic tail. (However in the geotail the particles 
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are also collisionless fce/orc the acceleration, which they are not in the corona). Consequently one 
cannot assign a temperature or uniform pressure to the plasma. Also, standard MHD resistivity 
becomes a questionable concept in these circumstances (even if one allows for “anomalous" 
resistivity). Lyons and Speiser ( 1985) have recognized this problem for the geotail and derived 
a formula for its “equivalent” resistivity, which can be negative in certain circumstances! It is 
a different concept from anomalous resistivity because collisions and wave-particle interactions 
play no role and consequently the magnetic energy is not converted into thermal particle energy 
(i.e. kinetic energy with no preference velocity direction) but into beam kinetic energy. 

A plasma under these circumstances is usually described by the Vlasov equation that deter- 
mines the evolution of the distribution function of the particles. However, in simulations for low 
density plasmas one often encounters difficulties here because densities tend to turn negative. 
Thus one is forced to model these collisionless current sheets at the most elementary level: by 
following each individual particle in its orbit, and considering the currents and magnetic field- 
generated by the these particle motions. Such simulations have been carried out for the geo- 
magnetic tail (see the review by Hockney and Eastwood (1988), Chap. 9-4). but not for the 
sheet in two-ribbon flares. The basic model and selfconsistency relations for these simulations 
will be discussed in Section 4. It will be shown that the numerical results for the geomagnetic 
tail can be scaled to those for the current sheet in a two-ribbon flare, since the dimensionless 
governing equations have identical form with similar numerical values for the constants. Hence 
the demonstration of the formation of a quasi-stationary current sheet and proton acceleration 
for the geotail (Hamilton 1976, 1980) is also valid for the corona. 

The final Section of this paper summarizes the main conclusions and suggests future obser- 
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vations. 


2 The Particle Acceleration Mechanism 

To explain the observed hard X-ray emission from flares as the result of the impact of a particle 
beam on the dense chromosphere, the mean energy of the particles has to be about 30 KeY if 
the beam is composed of electrons (Lin and Schwartz 19S7), and 200 KeV if protons are the 
carriers of the beam energy (Simnett 19S6). After acceleration in the current sheet, the mean 
free path (m.f.p.) for momentum loss by the electrons as well as the protons would be about 
lO 10 cm, equal to the length of the whole erupting structure. (\Ye have corrected here for the 
omission of a factor (ntt/nip) 1 / 2 in Eq. (2) of Martens (19SSa) for the m.f.p. of protons. Explicit 
expressions for the calculation of the m.f.p. ‘s are given in Trubnikov (1965). Other quantities 
discussed below are correctly calculated in Martens (19*Sa). and these calculations shall not be 
repeated here.) 

Given the particle energies, and an estimate of the electric field strength in the sheet of 2 
V cm -1 , found from the observed photospheric field strength and ribbon velocities (Kopp and 
Poletto 19S6). one can easily calculate the acceleration length of the particles in the sheet, by 
equating the required particle energy to the potential drop in the sheet over an acceleration 
length. The result for the solar sheet is an acceleration length of the order of 10 s cm, 5 orders of 
magnitude smaller than the m.f.p.! The actual m.f.p. of the particles may be 2 or at most 3 orders 
of magnitude smaller than 10 10 cm, because of collective (e.g. wave-particle) interactions, giving 
rise to anomalous resistivity, but an increase of 5 orders is unlikely for anomalous processes. 
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Clearly then, particles do not interact during their acceleration process, and one has to 
consider a different way to limit the acceleration path to about 10 5 cm. This situation is quite 
similar to that in the geomagnetic tail, where also the m.f.p.'s for protons and electrons are many 
orders of magnitude larger than the tail’s dimensions (see Lyons and Williams 1984 for a review). 
Early model calculations (Speiser 1965, 1967) prompted by in-situ satellite measurements of 
electric and magnetic fields have shown that in the geotail the orbits of individual particles are 
determined by the component of the magnetic field perpendicular to the sheet (B p ): after about 
half a gyration about this component the particles are ejected from the sheet with a very small 
pitch angle. Protons and electrons acquire the same velocity, and the beams emanating from 
the sheet are consequently neutral. For typical solar conditions this gyration radius is indeed 
about 10 5 cm for protons (Martens 198$a). 

Therefore the results of the numerical calculations of Speiser (1965), and in particular the 
approximate analytical expressions derived there, can be applied to the solar case. The results 
of this analysis (Martens 1988a) are: 

1. The drift of plasma into the sheet continuously supplies the needed 4 x 10 35 protons/sec. 

2. The sheet magnetic field (about 200 Gauss) generated by the proton current is acceptable 
despite the large number of accelerated protons. This stands in contrast with the result of Hoyng 
el al. ( 1976) for electron beams in single loops. 

3. Neutral particle beams are injected at the top of the postflare loops with almost zero pitch 
angle. 

4. The generation, timing, and spatial distribution of 7-rays and Hard X-rays, can be explained 
from a single power law proton spectrum, with cutoffs at about 200 KeV and 20 GeY. 



5. The effective resistivity in the current sheet, defined as the ratio of electric field to current 
density, is over 5 orders of magnitude larger than the classical Spitzer resistivity. 

Crucial in the derivation of the above results is the large electric field strength of 2 V cm" 1 . 5 
orders of magnitude larger than the Dreicer field. Kopp and Poletto derived this estimate from 
the l'xl field determined by the lateral diplacement of the flare ribbons across the normal 
component of the photospheric magnetic field. It was shown by Forbes and Priest (1984) that in 
a system with an ignorable coordinate along the filament axis (as is indicated by observations), 
the electric field at the x-type neutral point is exactly the same as near the ribbons for stationary 
reconnection. Note that this electric field is perpendicular to the magnetic field, except at the 
neutral point, where there is no magnetic field. 

This result is in contrast with the results of experiments in laboratory plasmas, where the 
electric field is usually limited to the Dreicer field (e.g. Stenzel el al. 1983. 1985). On the 
other hand, the electric field in the geomagnetic tail, determined by in situ observations of the 
reconnection rate, is about 8 orders of magnitude larger than the Dreicer field. In the experiment 
of Stenzel et al. the net electric field in the neutral sheet consists of two parts, an inductive 
part generated by the magnetic reconnection, and a space charge part, generated by charge 
build up at the anode and kathode at the endpoints of the sheet. The space charge electric field 
counteracts the inductive field in such a way that the net field is limited to the Dreicer field. 

The current sheet in a two-ribbon flare and the geomagnetic tail are not bounded in the 
way a laboratory plasma is, and therefore charge separation at the endpoints of the sheet is not 
possible. For example, charge separation at the endpoints of the flare current sheet (see Fig. 1) 
would immediately be cancelled by removal of the excess charge carriers along the fieldlines of 



ORIGINAL PAGE IS 
OF POOR QUALITY 

the postflare loops connecting to the for all practical purposes infinite reservoir of particles in 
the chromosphere and photosphere. Consequently the net electric field is completely determined 
by the reconnection rate of the magnetic field. 

As was noted by Kopp and Poletto (1986) the strict mathematical result of Forbes and 
Priest (19S4) for the inductive electric field is not applicable in the absence of perfect symmetry 
along the axis of the filament. However, according to Faraday’s law a net electric field is always 
induced along the magnetic null line when two flux systems exchange magnetic flux, irrespective 
of symmetry. Farady’s law can be integrated along the x-tvpe neutral line, as done in Kopp and 
Poletto (their Eq. (4)) to express the net voltage drop along the neutral line as a function of the 
increase in net flux in the postflare loop system. The latter is equal to the photospheric line of 
sight flux bypassed by the lateral motion of the flare ribbons, that is directly measured. The two 
flux systems in the present model are (1) the loops overlying the filament, and (2) the postdate 
loops plus the field of the filament (see Fig. 1). The flux exchange leads to the formation of n<*w 
postflare loops through reconnection within the current sheet. 

Clearly then, large electric fields in the current sheet that is created in the wake of an 
erupting filament can only be avoided when either there is a space charge electric field, or there 
is no reconnection at the top of the postflare loops (no flux exchange). The first possibility wa> 
ruled out above. The second one would imply that the same arcade of loops connects to the 
x-tvpe neutral line during the flare, and its footpoints would not change position, in contrast to 
the observed lateral motion of the ribbons. However, it is possible that the observed Ho ribbons 
have nothing to do with loops connecting to the current sheet. Instead some MHD instability 
can be envisaged to take place at sequentially higher and higher loops that already existed before 
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the flare. The ribbons then still are the chromospheric signature of the energy release within 
each subsequent arcade of loops. In this picture the whole preflare structure would simply be 
stretched out during the flare, without any change in topology. 

There are several possible observations that could determine whether the reconnection or 
the stretching picture is correct. The first one is by directly measuring the electric field during 
flares: a field of several V cm' 1 would imply reconnection. A second straightforward test is 
the observation of changes in topology. Suppose one can observe a preflare loop overlying the 
filament before eruption, and determine its footpoint positions. If, during or after the eruption, 
the ribbons connecting the postflare loops cover or go beyond the footpoints of the preflare loop, 
reconnection must have taken place. 

We are not aware of any direct observation of topology changes, but electric field measure- 
ments have recently been attempted. Foukal et al. (1983. 1986. 1987) have derived an upper 
limit of about 10 Y cm' 1 for the macroscopic electric field in eruptive prominences and post flare 
loops by observations of Stark broadening in Balmer and Paschen lines. Future observation?- of 
the polarization of the Paschen lines may allow the measurement of macroscopic electric fiekl> 
of substantially less than 10 V cm -3 , and thus the results of Kopp and Poletto may be verified. 

We conclude that the presence of super Dreicer fields in two-ribbon flares is not unambigu- 
ously demonstrated from observations, but does provide the most straightforward interpretation 
of the observations of the formation of postflare loops, and is also suggested by the presence of 
super Dreicer fields in the similar circumstances of the geomagnetic tail. 

A three way comparison between the physical situation in the current sheet of a two-ribbon 
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flare, the geomagnetic tail and the laboratory plasma experiment of Stenzel et al. ( 1983 , 1985 ) 
is presented in Table 1. The first part of the Table gives the observed quantities and their 
references. From these, and assuming that the beam generation mechanism of Speiser ( 1965 . 
1967) is at work, one can calculate the sheet physical parameters such as the drift velocity, the 
m.f.p. ’s, the proton acceleration length, the effective resistivity, number of accelerated protons, 
net current, etc. (see Martens 1988a for the results for the geotail and the corona). 


The effective resistivity is a parameter of particular significance since it is more than 5 orders 
of magnitude larger than the classical Spitzer resistivity in the coronal current sheet. It is derived 
by considering the momentum equation for the protons in the current sheet, 


c/1 - 7?i r l -1 

m r — = e£ — - — . 

r dt m .f.p. 


( 1 ) 


In a stationary situation the left hand side of this equation vanishes on the average. The proton 
acceleration inside the sheet is limited by whatever mechanism ejects the protons from the sheet. 
Using the definition of the proton current density, j = n p e\' r one finds for the effective resistivity 
defined above 


’hJJ = 


n 'r ] p 


v p e 2 m.f.p.' 


(2j 


The m.f.p. for proton-proton collisions is very large, so collisions are a very ineffective limiter. 
Speiser’s (1965) calculations show that the gyration of the protons along the small remaining 
perpendicular component of the magnetic field results in the ejection of the accelerated protons 
after half a gyration. This limits the acceleration path of the protons to roughly 2 gyroradii of 
B p , much less than the collisional m.f.p. Replacing the m.f.p. by two gyroradii one finds 


’hff = 


B p 

n p ec' 


( 3 ) 
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A more detailed derivation of and in particular its dependence upon the initial values of 
V p is given in Lyons and Speiser (1985). A peculiar aspect of their results is that protons with 
very high intial velocities are actually decelerated in the sheet, and hence their 7? e // is negative. 

More recently (Martens and Young 1988) we have derived some new results for the current 
sheet in two-ribbon flares that strengthen the conclusions above. The thickness of the current 
sheet (£) may be estimated from the induction law with the assumption of a quasistationary 
sheet, which is based on the observation that the duration of a flare is much longer than the 
Alfven travel time through the sheet. The result is 

TjC ^ 

f — - — — ft 430cm. 1 4 ) 

2-\ d 

where Yj is the plasma drift velocity into the sheet (ft 15 km sec -1 ). The net energy conversion 

I 

in the sheet (P) can be derived as follows. 

P = ILL ft 6 x 10 2S er/y/#ec. (5) 

where I is the current in the sheet. E the electric field, and L the length of the sheet. This result 

is in good agreement with the energy requirements in the impulsive phase in large flare-. 

i 

Finally, from magnetic flux conservation we derive an estimate for the perpendicular magnetic 
field component. Equating the inflow of magnetic flux to the outflow, one finds 

^ in B sheet = Vovt^p- ( 0 ) 

The inflow velocity is equal to the drift velocity, but the outflow velocity is not known a priori. 
Its maximum value is the local Alfven speed. For the geomagnetic tail all the quantities in Eq. 
( (5), except have been measured (Lyons and Williams 1984), and 1 ’, ut is found to be equal 
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to the Alfven speed. This corresponds to maximum energy conversion in the sheet. We assume 
that the same situation holds for the flare current sheet, and the result is B p ~ 0.7 Gauss. Given 
this value the average energy of the accelerated protons can be calculated with the approximate 
formula given by Speiser (1965), 

Up « 2 m v c 2 (E / Bp) 7 ss 174A'eV. (7) 

This result is consistent with the average energy of 200 KeV required from observations (Simnett 
19SG). that was used as an input parameter in our calculation.-.. Hence the assumption V = 
1^1 is justified by the observational constraints on our model. We further note that 1 '• .. r = ]' A i.» 
a common characteristic of many different sorts of sheet models (see Priest 19S 1 . Chap. 3. for 
a review). Note also that the B p calculated is the one at the edge of the sheet. Near the center 
B v is much smaller, even reversing sign at the axis. Consequently the protons accelerated near 
the center of the sheet can acquire much higher energies, theoretically up to the full potential 
drop along the sheet of 20 GeY. 

In the second part of Table 1 the results for some dimensionless parameters of the three 
plasmas are compared. It is obvious from the Debye number, plasma- J. and ratios of plasma- to 
gyrofrequency that the plasma of the geotail and the flare current sheet are quite similar, despite 
their greatly different dimensional parameters. Equally obvious is the great difference with 
StenzeTs plasma experiment, despite the greater similarity in dimensional parameters between 
the coronal current sheet and the experiment. The sheet thickness for the geotail and coronal 
current sheet are calculated with Eq.( 4) above. For the geotail the result is in agreement with 
spacecraft measurements, while in the corona direct experimental verification is impossible at 
this time. The thickness of the laboratory sheet is directly observed. In StenzeTs laboratory 
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current sheet most of the current is not carried by protons, but by electrons. The net electric 
field is smaller than the Dreicer field because of the charge build-up at the anode and kathode, 
that cancels the inductive electric field. The derivation of the perpendicular field component for 
the coronal current sheet (Eq. ( 6)). only indirectly verifiable, was discussed above. The effective 
resistivity (Eq. ( 3)) can be verified indirectly from the observed net energy conversion rates 
(Eq. ( 5)), giving good agreement for the geotail and the corona. For the laboratory experiment 
rj e jj is found as the ratio between current density and electric field. 

Our conclusion from the comparison in Table 1 is that the particle acceleration process in 
the geomagnetic tail and in the current sheet of a two-ribbon flare, both astrophvsical plasmas 
without walls and with resembling dimensionless characteristics, are likely to be similar. 


3 Observational Signatures of Neutral Beams 

Protons have been proposed as the primary means of energy transfer in flares from early on 
(Carmichael 1964. Elliot 1964. Schatzman 1967). but in the last decade or two electron beams 
have been the main focus of attention. Recently Simnett and Strong (1984) and Simnett (1986.1 
have argued that the electron beam hypothesis should be reexamined, because it fails to explain 
some key observations, and that neutral beams with protons as the carriers of the beam energy 
provide a straightforward explanation for the same observations. Their main arguments are: 
(1) SMM-FCS observations of the Ca XIX ion (Antonucci et al. 1982, 1984) show sometimes 
large non-thermal line broadening and high upflow velocities prior to the onset of the impulsive 
phase (hatd X-ray emission) of flares. This requires large energy deposition in the chromosphere 



prior to the impulsive phase, and electron beams cannot do this without generating substantial 
thick target X-ray emission, which is bv definition not observed. Proton beams can generate 
such heating without producing X-rays. Only when the heating rate in the chromosphere due to 
the bombardment by protons surpasses a certain treshold value, a superheated plasma will be 
generated (Fisher et al. 1985, Fisher 1987) that emits hard X-rays: this is the impulsive phase 
of the flare. (Note that in this scenario the X-ray emission is generated by thermal electrons). 

(2) X-rays generated by thick target impact of electron beams should be anisotropical and 
polarized, in contrast to thermal X-rays. Observations show isotropy (Kane et al. 1980) from 
simultanuous stereoscopic data, and only small upper limits to polarization (Tramiel et al 
1984) incompatible with electron beams. The data are fully compatible with the generation of 
a superhot plasma in the chromosphere by proton beams. 

(3) Observations with the SMM-GRS instrument (Chupp et al. 1982. Forrest ct al. 19S1) 
demonstrate conclusively that high energy protons (> 30 MeV) must be present from the onset 
of the impulsive phase of “ -ray flares: the * -rays are cotemporal with the hard X-rays within the 
time resolution of the GRS instrument (2 sec.). The number of high energy protons is at least 
two orders of magnitude larger than that of high energy electrons (Ramathv and Murphy 1984 j. 
One acceleration mechanism for both MeY and KeY protons (as proposed in this paper) seems 
more straightforward than two different but exactly simultanuous mechanisms for KeY-electrons 
and MeY-protons. 

Valid proof for the existence of proton beams, however, can never be obtained from arguments 
against electron beams or plausibility arguments alone. Dennis (1988) has compared proton 
beams to the Loch Ness monster just because all the evidence for their existence is of the above 
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nature. The observational difficulty with hecta KeV protons is that they do not generate X-rays 
upon impact, or n - waves during their transport, as deka KeV electrons do. 

We will discuss three independent observational tests for proton beams here, two giving 
evidence for proton beams, and the other inconclusive. 

Orral and Zirker (1976) and Canfield and Chang (1985) have investigated the interaction 
of proton beams with the parly ionized chromosphere. They find that the main source of La 
wing emission is charge exchange from the proton beam followed by spontaneous bound-bound 
transitions of the beam atoms. The increase in La wing emission that is predicted for typical 
beam fluxes and proton energies is several orders of magnitude larger than what has been 
observed from Skylab for two flares (Canfield and VanHoosier 1980). but ‘it is not at all clear 
that those observations were obtained at a time when we expect protons to be present during 
these events" (Canfield 19S8). Hence no clear indication or counterindication for proton beams is 
found here. Future space based measurements with UV spectrometers with a spectral resolution 
of 1 Vi or less can resolve this question 

Henoux and C’hambe (19SS) calculated the polarization in the line emission of hydrouen 
atoms •excited by a beam of protons or electrons. They found that the direction of vibration of 
the electric field vector generating the the Ho emission is polarized parallel to the incident beam 
for dekaeV electrons and deka KeV protons. (The polarization depends mainly on the vdocily 
of the beam panicles, hence the m p /m e energy ratio). Assuming beam incidence perpendicular 
to the solar surface, an earth bound polarimeter would observe maximum intensity in the flare 
site - center of disk direction. For deka KeV electrons and deka MeV protons the electric field 
vector is polarized in the plane perpendicular to the beam, and hence an earthbound observer 
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sees polarization perpendicular to the flare site - disk center direction. 

Three chromospheric flares were observed at Meudon on July 11 1982 (sec the figures in 
Henoux 1989, these proceedings), and the results clearly show polarization in the disk center 
direction for several Ha emitting spots. The timing coincides with the period from the onset to 
the peak of the soft X-ray emission as observed by GOES 2, which corresponds to the impulsive 
phase. Henoux et al. (1988) rule out deka eY electrons for two reasons: (1) the m.f.p. of these 
electrons is so small that they have to come from the transition region, which is unlikely, and 
more important (2) these electrons do not carry enough energy to power the simultanuou- X- 
rays. Deka KeV electrons, and hence standard electron beams, are ruled out because they would 
cause polarization in the perpendicular direction to the one observed. 

For typical preflare coronal values of T & 2 x 10 G I\ and n t % 2 x 10 9 cm -? '. and the height oi 
the current sheet at the onset of the flare of about 10 9 cm (Martens and Kuin 19S9. Takakura < 1 
al. 19S6).one finds a prefiare coronal column density of about 4 x 10~ € ' g cm~~ . less than model 
F] in Henoux et al. (19S8). From extrapolation of their results we find thai the energy flux 
in the neutral beam that is required is to produce the polarization is about 10 10 e rg ar<~- »“•. 
and the protons need to have an initial energy in excess of 200 KeY to reach the Ho producing 
layer. This is in good agreement with the theoretical results of the previous section: there i> a 
lower energy cutoff in the energy spectrum at 200 KeY. and with the power calculated in Eq. 
( 5), and a surface area of the ribbons of about 10 19 cm 2 (length 10 10 cm. width 5 x 10 s cm ). the 
predicted energy flux in the beams is 10 10 erg cm -2 s~ l 

Similar observations of impact polarization in the S I 1437 A line have been made with the 
SMM-UVSP instrument (Henoux tl al. 1983). They found polarization up to 25/7 within 3<letr 
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of the flare site - center disk direction, again implying impact by electrons in the deka eV range 
or protons in the deka KeV range, and again totally inconsistent with the polarization direction 
expected from deka KeV electron beams. A polarization of 9% was observed during the peak 
of the impulsive phase (the onset of the soft X-rays), and 25% polarization occurred during the 
decay phase of the flare, with no polarization in between. While the latter peak in polarization 
is most easily accounted for by impact of deka eV electrons carrying down heat flux from the hot 
and dense postflare loops as was proposed by Henoux ei al. (1983). the polarization at the onset 
of the flare cannot be explained this way. in the absence of preflare dense loops. Neutral beam? 
on the other hand, penetrate easily through the tenuous preflare plasma, as was demonstrated 
above. We thus submit that the flare on?et polarization is caused by neutral beams, and that 
this effect vanishes as the flare loops become denser, because the proton m.f.p. in the corona 
declines. Later on, in the postflare phase, heat flux carrying electrons revive the polarization. 
Clearly more observations are required to test this scenario. 

Finally we will briefly discuss an old observation that supports the specific proton acceleration 
mechanism described in the previous section. The maximum attainable proton energy is defined 
by the full potential drop along the current sheet. With the electric field strength of 2 Y cm -1 
and sheet length of about 10 10 cm this maximum energy is about 10 GeY. This should be 
regarded as an absolute upper limit, since only the largest two-ribbon flares have such long 
filaments. Observations of high energy solar flare protons with the worldwide neutron monitor 
network in the early 70’s (Heristchi el al. 1971. 1976). have produced direct evidence for a high 
energy cutoff in the proton spectrum in the range from 2 to 5 GeV for seven flares for which the 
cutoff energy could be calculated, close to our theoretical value. 
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4 Numerical Models 


4.1 Single Particle Orbits 

We have calculated the trajectories of protons through a model two-ribbon flare current sheet. 


The equation of motion for a collisionless proton is 


dV ,r. VxB , 

m P — = e £+ ). 

h dt c 


We use a coordinate system with y along the axis of the sheet, i in the direction of the width 
of the sheet, and x along the thickness. Consequently the variation of the magnetic field in 
the i direction is strongest. We assume a stationary sheet (d t = 0) with translation symmetry 
along its axis (d y = 0). as in Speiser (1965). The induction electric field is directed along y and 
constant in space and time, by virtue of these assumptions, the absence of charge separation, 
and Faraday ? s law. From Table 1. E y = 2 V era" 1 . The i component of the magnetic field i- 
constant outside the sheet and varies linearly inside the sheet. 


B : = B s for r < -(• 

B : = f or x < f 

c 

B : = -B s for x >6. (9) 

with B a = 200 Gauss (Table 1) and 6 = 430 cm (Eq. ( 4). The y component of the magnetic field 
is zero everywhere, and the x component, constant in space and time, is found from Eq. ( 6). 
B p = 0.7 Gauss. In reality B p changes slowly in the i direction, reversing sign at the sheet axis. 
Here we use the B p at the edge of the sheet. The numerical results show that the excursion of 
a proton in the f direction inside the sheet is about equal to the gyroradius calculated from B,. 


19 


and the final velocity of the particle, as 10 5 cm. This is small compared to the scale for variation 
in B p (10 9 cm), and hence the constancy of B v is a good approximation. 

Outside the sheet the magnetic field gradient vanishes, and the particle motion is fully 
Integra ble (Sivukhin 1965). The particle motion then consists of three components: (1 ) constant 
velocity along the magnetic field, (2) gyration around the magnetic field lines. (3 ) ExB drift into 
the sheet. Our numerical results reproduce this correctly. The interesting part of the trajectories 
is inside the sheet, where the equation of motion is not fully integrable. For most calculations 
we have restricted ourselves therefore to the time between entry and exit of an orbit. Figure 2 
displays the Tesult of one of our calculations. 

Inside the sheet the length scale for the variation of the magnetic field is comparable to the 
proton gyroradius. and the three components of the particle motion described above couple in 
a nonlinear way. We found that the protons indeed follow trajectories of the type described by 
Speiser (1965. 1967. 1968). as was surmised in Sect. 2: the protons oscillate rapidly accros t he 
midplane of the sheet, at the same time are accelerated along the electric field, and finally are 
slowly deflected by the small perpendicular component of the magnetic field. After half such a 
gyration they are ejected from the sheet in a direction almost parallel to the magnetic field. For 
protons with initial energy of the order of the typical coronal thermal energy, the approximation 
formula for the final particle energy (Eq. ( 7)) was found to be correct within 2% for each of 
many test orbits. 

Despite the nonlinearity of the equation of motion some analytical results can be established 
that lead to insight in the numerical results, and that are useful in the derivation of a fully 
selfconsistent model. In fact. Speiser (1965) derived an analytical approximation to the proton 
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trajectories inside the sheet. 


Because of the absence of y dependence, the magnetic vector potential can be reduced to 
one component, 

J9 = V x (A y y). (10) 

In the particular case of the model calculation of this section 

B x 2 

. * Ay — — Bp z — 2 By t- (11) 

The y component of the equation of motion (Eq. ( 8)) may be integrated to yield the conserved 
canonical momentum in the y direction. 

P y = m p Y y - C -^ = Cj (12, 

The constancy of B v yields yet another integral 

m V. - = C 2 . (13) 

y " C 

Finally, the stationary gives rise to the third integral, expressing conservation of net energy. 
Internal multiplication of Eq. ( 8) with V and integration yields 

- eE y y = C 3 . (in 

The three equations above describe the velocity vector field completely as a function of posi- 
tion and initial conditions. However, the constants of motion are not in involution (Chen and 
Pahnadesso 1985) and an analytical solution cannot be derived. The first and the last equation 
are valid irrespective of the form of A y , the second depends on the constancy of B p . However, 
it should be noted that for practical purposes B p is usually constant for a typical trajectory. 
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For constant B p one can further simplify these equations by changing to a frame in which the 
electric field is transformed away (see Speiser 1968 for further details). 


One cannot derive the ejection energy and ejection angle from the above conserved quantities 
alone; the numerical integration is needed. However, if we use the numerical result that the 
particles with low initial velocities are ejected almost parallel to the magnetic field, so that 
|Vy| k |V/ t ,|, one can derive from Eq. ( 13) and Eq. ( 14) by eliminating y. 


m-1 7 


2 m p c\fY, 


(IV. 


confirming the estimate of Eq. ( 7). The qualification low initial velocity refers here to the fin;, 
velocity V/ = 


4.2 Selfconsistent Sheet Models 

In a selfconsistent model one has to take into account that the currents are the result of the 
relative motion of the protons and the electrons, in other words, the net currents that are found 
by considering the summation over the individual particle orbits calculated in the previoii- 
section need to be consistent with the currents in the originally assumed magnetic field. (In the 
model of Sect. 4.1 j y = |!rf inside the sheet, and the current density vanishes outside the sheet.) 

It may be shown that in a current sheet of the type of the geomagnetic tail the curreni is 
mainly carried by the protons, because of their larger gyroradius. and that the electrons will 
move in such a way as to preserve charge neutrality, and prevent currents other than parallel to 
the sheet (Eastwood 1972, 1974). In this case the basic equations governing the formation and 
evolution of a sheet are: 
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(1) Faraday's law 


1 d§ 

c dt 


= -V x E, 


(2) Ampere's law 

V x B = — J. 

c 

(3) the momentum equation for the protons. 


dV p VxB x 

" , '>rfT = <l£+_ r ) ' 


(4) the definition of the current density 


j ■ E — n. f. V • E. 

j x £ = (7. 


( 10 ) 


(IT) 


( 1 *.. 


( 10 ) 


where the bar indicates the local average over the velocity of all protons at that location, u 
is the proton density, found from summation over the individual protons. Note that since this 
density is found from summation mass continuity is automatically satisfied, and there is ru> 
risk of the density becoming negative at some point in space, as often happens in numerical 
simulations for low density plasmas with the MHD equations (see Hockney and Eastwood 19S y . 
Chap. 9-4 for a detailed discussion of this approach). 

We now consider a sheet of the same type as in Sect. 4.1, i.e. the same coordinate system. 
B y = 0 and d y = 0. We drop for the moment the stationaritv, and allow for variations in the 
i direction. Using the above constraints and Faraday’s law, Eq. ( 17). it is found that E in 
the x,z-plane is a gradient type vector. However, because of the absence of charge separation 
(V • E = 0) and vanishing E at large distances, it can be shown that E x = E : = 0. and hence 
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Eastwood’s (1972, 1974) result that j T = j s = 0 is rather plausible. Thus we have. 

J = (0.j*,0) 

E = (0,£ y .0) 

A = (0,.4„.0). (20) 

Note that these quantities are not independent of time and the r coordinate, as in the previous 


section. Ampere’s law now reduces to 


V : Ay = --Jy. 


where V 2 = # 2 4- d: in the absence of y-dependence. 


Of the constants of motion for individual particles only one remains, the canonical momentum 


in the v-direction 


, , . e .4 „ 

= m p \ y — = con-tinni. 


Eqs. ( 19). ( 21), and ( 22) may he combined to give the field equation for the problem 




v — 

c- / <- e 


with iJp r denoting the (variable) proton plasmafrequency. 


4tttj r c 


The field equation has to be solved in conjunction with the two remaining components of the 
momentum equation for the individual proton orbits, which after the elimination of V y can be 
written in the following convenient notation 


dV ^1 j/e.4 v \ 2 eA 
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Note that A y is a function of time in this equation, so energy is not a conserved quantity here. 

The above equations model the three dimensional motions of the protons, laking into account 
the effect of finite gyroradius and frequency, that are absent in the MHD description. The 
model also represents sound- and Alfven waves. Numerical simulations have been done for the 
geomagnetic tail (Hamilton 1976. 1980. see Hockney and Eastwood 1988. Chap. 9-4, for an 
overview) and they demonstrate clearly the formation of a current sheet when particles arc 
injected into an originally current-free x-type magnetic field configuration. It is further found 
that the sheet configuration is unstable for a cross-sheet field component (D p ) smaller than a 
certain threshold. It may be expected that the current sheet in a two-ribbon flare is unstable 
for this reason, because of its large D s to B P ratio in comparison to the geotail (Table 1). 

Wo note that there is a natural length scale in the field equation ( 23). ( = that will 
probably determine the thickness of the current she*et. This length scale is only dependent on 
the density through Eq. ( 24). With the typical values of Table 1 we find ( * 10' cm for the 
geotail, and ( k 300 cm for the solar current sheet, in both cases close to the thickness of the 
current sheet, obtained from the induction equation ( 4). Comparison between the expressions 
for ft and i shows this equality is not at all trivial. This fortunate coincidence further suggests 
that the results for the evolution of the current sheet in the geomagnetic 1 ail can be scaled 
directly to the current sheet of a two-ribbon flare. 

We have formalized this by considering the dimensionless model equations. First we intro- 
duce the following parameter and units: 

. _ El 

B s ' 
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to 


m p c 
eB s ' 


V D = 


E o£ 

B„ ‘ 


(2G) 


Hence the time unit is the inverse proton gyrofrequency, the velocity is normalized such that the 
final velocity is of the order unity, and the initial (thermal) velocity and drift velocity are very 
small for both the geotail and the solar corona. The parameter e expresses the dominance of 
the sheet magnetic field over the the perpendicular field. The dimensionless equations of motion 
become 


y = sh(x. z.1)+ ezg{r.:.1) - 

~ = (27) 

The functions f.g-h represent respectively B : . B~. and E v . and form the partial derivatives of 
the dimensionless potential .4(.r. r. /). By construction they are all of order unity. The s values 
for the geotail and the corona are given in Table 1. 


Having introduced a unit for time and velocity, the unit for length is automatically deter- 


mined 


tn F c 2 £ 0 

Lo _ \. oto _ . 


(2s ) 


Introducing tills in the field equation ( 23), and using the obvious definition 


P — 

■* Vi° ~ c ’ 


(29) 


one finds 


(-V 2 + an)A = anP 


(30) 
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Here a new parameter o is introduced. 


q = 


(B p B s )i 


For the geotail the value is 0.33, for the two- ribbon flare 0.43. 


(31 ) 


Thus the only two relevant parameters in this problem, e and a, are quite similar in value 
for the geotail and the corona. Hence the numerical results for the evolution of the geomagnetic 
tail can simply be scaled to those for the corona, with the use of the physical units given above. 
The results of Hamilton (19S0) show that for an c-value as low as that calculated above for the 
coronal sheet the reconnection process will be non-stationary. with intermittent bursts, and the 
generation of multiple X- and O-type neutral lines in the sheet. This is not inconsistent with 
the observed nature of the hard X-ray emission. 


5 Conclusions 

We have argued that there is a great similarity between the physical situation of the geomagnet ic 
tail and the current sheet in the wake of an erupting filament. We have demonstrated that t h** 
dimensionless parameters characterizing these plasmas are comparable, and that the prediction 
for total energy conversion, number of accelerated particles, and proton energies for a two-ribbon 
flare are in good agreement with observations. 

Also, from a numerical computation of individual proton orbits we have verified the applica- 
bility of Speiser’s (1965. 1967) proton acceleration mechanism to a sheet with solar parameters. 
Finally we have shown that the similarity between the solar current sheet and the geotail goes 
so far that the dimensionless evolution equations for both systems, under the assumption of 
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invariance along the sheet axis, and dominance of proton currents, have similar constants. Con- 
sequently the results of numerical simulations of the geotail (Hamilton 1980) can be scaled to 
the corona. 

Further we have reviewed the observational material pertinent to neutral beams populated 
with hekta KeV protons, and standard deka Kev electron beams. We find that in addition to 
Simnett's (1986) arguments for neutral beams, recent line polarization measurements (Henoux 
el al. 19S9. these proceedings) seem to rule out standard electron beams in favour of neutral 
beams. However, the observations are too few. and their interpretation is not yet sufficiently 
established. 

Future measurements ofline polarization - either ground based Ha or W space observations - 
are clearly desirable, as well as high spectral and temporal resolution measurements ol red shifted 
Lymann. advocated by Canfield and Chang ( 19S5). On the theoretical side, the application o! 
codes developed for the geotail to solar type sheets is called for. 

Studies of coronal plasma processes at the Harvard-Smit hsonian Center for Astrophysics an 
supported by .Y.4S.-1 grant XAGW'-llH. 
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Tables 
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Table 1: A Comparison between the Plasmaphysics of the Geomagnetic Tail, the Current Sheet 
in a Two-Ribbon Flare, and Stenzel’s Current Sheet Experiment. 


|| Physical Quantity Geomagnetic Tail | Two-Ribbon Flare Laboratory 

|| Observations Hi 

Sheet Magnetic field 

2 x 10-“ Gauss (1,2) 

200 Gauss (3) 

(typical for active region) 

12 Gauss (14) 

Sheet Electrical Field 

2.5 x 10 _e Volt/cm (1) 

2 Volt/cm (4,5) 

0.1 Volt/cm (14.15) 

Energy of Accelerated 
Protons 

2-100 KeV (1,6) 
(power law) 

200 Kev - MeV range (7) 
(inferred from Hard X-Rays) 

N.A. 

Sheet Length 

2.6 x 10 1U cm (8) 

1.2 x 10 1U cm (9) 

(inferred from filament length) 

180 cm (14) 

Sheet Width 
thickness!) 

1.3-13 x 10 IU cm (8) 

10* cm (9,10) 

(inferred from width ribbons) 


Density of Ambient 
Plasma 

0.02-5 cm _iJ (11) “ 

10 10 cm- 15 (12) 

(typical for active region) 


Temperature of 
Ambient Plasma 

1.5 x 10* K (13) 

2.0- 3.0 x 10* K (12) ““ 

(typical for active region) 

m 

Hiffittuiun 

|( Dimensionless Quantities 

Number of Particles 
in Debye Sphere 

4.0 x 10 1J 

3.9 x 10 7 

3.4 x 10“ 

Plasma Beta (outside sheet) 

5.2 x 10" J 

4.3 x 10"^ 

2 

Plasma Freq. /Gyrofreq. 

9 

1.6 

267 

Sheet Thickness/Proton 
Gyroradius before Acceleration 

34(1) 

52 

0.83 (14) 

Sheet Thickness/Proton 
Gyroradius after Acceleration 

0.88 (1) 

1.9 

N.A. 

Electric Field/Dreicer Field 

1.2 x 10 s 

6.7 x 10* 

0.15 

Parallel /Perpendicular 
Compononent Magnetic Field 

20(1) 

290 

not measured 

Effective Resistivity/ 
Spitzer Resistivity 

1.4 x 10 9 

3.0 x 10 s 

25 

Observed/Calculated Proton 
Energy after Acceleration 

1.5 

1.2 

N.A. 
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Figure Captions 


Figure 1: The magnetic field structure of an erupting filament, and the associated two- ribbon 
flare, viewed in a cross-section perpendicular to the filament axis. The location of energy con- 
version processes is indicated. 

Figure 2: A Typical Proton Orbit in the Current Sheet of a Two- Ribbon Flare. Inside the 
sheet the motion along the y and f axis is scaled down by a factor 250 resp. 500 for clarity of 


presentation. 
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Scaling Laws 


o 


Summary 

Scaling laws are derived for the one-dimensional time-dependent Euler equa- 
tions that describe the evolution of a spherically symmetric stellar atmosphere. 
With these scaling laws the results of the time-dependent calculations by Ko- 
revaar (1989) obtained for one star are applicable over the whole Herzsprung- 
Russel diagram and even to elliptic galaxies. The scaling is exact for stars with 
the same M/R - ratio and a good approximation for stars with a different M/R- 
ratio. The global relaxation oscillation found by Korevaar ( 1989 ) is scaled to 
main sequence stars, a solar coronal hole, cool giants and elliptic galaxies. 

key words: galaxies: coronae of — galaxies: elliptical — hydrodynamics — 
numerical methods — stars: coronae of — stars: mass loss 



Scaling Laws 


3 


1 Introduction 

If a stellar corona is heated by a mechanical flux greater than a certain limit, 
it cannot exist as a stationary extended corona. This was found independently 
by Hearn and Vardavas (1981), Hammer (1982) and Souffrin (1982). Hearn et 
al. (1983) suggested that a corona heated by a mechanical flux above that limit 
will undergo a global relaxation oscillation. Time-dependent calculations by 
Korevaar (1989) have confirmed this. For the one global relaxation oscillation 
that has been calculated, the corona collapses and rebuilds in about 2 years. 
It is then quasi-stationary for about 5 years, after which a new collapse occurs. 
Korevaar (1989) give a detailed description of the relaxation oscillation and 
explain the physical time scales that dominate the different phases of the global 
relaxation oscillation. 

The corona collapses because thermal conduction, which is independent of the 
density, is unable to heat the outer layers of the corona sufficiently, when their 
density surpasses a certain limit. A density pulse is formed at about 50 stellar 
radii which grows in amplitude as it falls inwards. It hits the photosphere 
hypersonically with velocities of about Mach 25. Most of the momentum of 
this infalling pulse is transferred inelastically to the photosphere and much of 
the kinetic energy is radiated away. The rest of the momentum is reflected 
in a pulse with an amplitude 2 orders of magnitude lower than the infalling 
pulse. The rest of the kinetic energy of the infalling density pulse is converted 
to kinetic energy of the reflected pulse and heats the gas. When sufficient 
material has moved into the photosphere a new corona is formed by the me- 
chanical heating, and it expands rapidly outwards. Dtiring this expansion the 
photosphere oscillates temporarily at the Brunt- Vaisala resonance frequency, 
responding to the momentum transferred by the infalling density pulse. Waves 
originate from the photosphere and enter the corona. This temporary oscil- 
lation is damped by radiation and conduction. The inner corona stabilizes 
and the remaining waves travel outward until they are damped. The corona 
then becomes quasi-stationary until a new density pulse is formed in the outer 
layers and the relaxation oscillation reoccurs. 
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Korevaar (1989) calculated the relaxation oscillation for the hot supergiant 
C Ori with an effective temperature of T t g = 31000 A', a radius of R „ = 27.8 R © 
and a mass of Af» = 44.7 A/©, reduced to 25.1 Af© to include the effective reduc- 
tion in gravity through electron scattering. One cycle of the global relaxation 
oscillation requires about 2 hours on a Cyber 205 supercomputer. This makes 
it impractical to repeat the calculations for a whole range of stars. 

Fortunately this is not necessary. The one-dimensional Euler equations that 
describe the time-dependent behaviour of a spherically symmetric stellar at- 
mosphere can be made dimensionless. In the present paper scaling laws are 
derived that extend the applicability of the calculations of the global relax- 
ation oscillation over the whole Hertzsprung-Russel diagram and to elliptic 
galaxies. The scaling is exact between stars with an equal M/R - ratio and a 
good approximation for stars with a different M/R- ratio. 

Hammer (1984) has derived a similar set of scaling laws for the one-dimensional 
time-independent Euler equations. In the present paper the time dependency 
is included as well. Hammer’s (1984) scaling laws stem from a homologous 
transformation of the equations, while the scaling laws in the present work 
are derived by making the flow equations dimensionless. The two methods 
are essentially the same. Hammer (1984) replaced each physical variable / by 
/' • kj and then derived relations between the k/s for the various variables. 
Initially Hammer allows each kj to depend on /, but later he requires for a 
proper scaling that all kj's are constant. In the present paper this is required 
from the beginning. 


2 Scaling laws 

A set of equations describing the evolution of a given system is expressed in a 
consistent set of units, usually MKS or Gaussian in astrophysics. The actual 
choice of units is arbitrary, and one may often employ this freedom to reduce 
considerably the number of physical parameters in the system. For example, 
by choosing the unit of mass in stellar evolution calculations as A/„ in stead 



Scaling Laws 


•5 


of kilograms or grams, the system of equations becomes independent of this 
parameter. By reducing the number of parameters one also reduces the number 
of numerical calculations required to investigate the solutions over the range 
of the parameters. 


In most cases it is impossible to eliminate all parameters from the equations. 
The ones that remain lead to scaling laws, so that for a constant value of this 
parameter one finds identical solutions. The relation between the parameter 
and the original stellar parameters is the actual scaling law. Two examples 
are seen in what follows, the expressions for the parameters / and a. 


The Euler equations for a one-dimensional time-dependent stellar atmosphere 
with spherical symmetry are 


1 


V r 


dt} pTVr) + X* dx T 

1 d 


e '=2 * r ' + ->-r 

TlTr 

Pr — Pt i 

P 

d 1 9 , 2 \ . 

K* + 7}K = °- 

i d 


( 2/ . \\ GM* %Pt 

(X r (p r t-V+Pr)j =-Pr-^j- + — • 


a io/ 2 , s\ GM* ~ 

+ + P'*) = -P'” ~ + Q 

Qtoi _ Qcond + Qmtch + Qtad + Qop 


tot 
r * 


(«) 

(b) 

(c) 

id) 

(e) 

if) 


A detailed description of these equations has been given by Korevaar and Van 
Leer (1988). The dimensionless Euler equations are derived below. Dimen- 
sionless quantities are indicated by symbols without a subscript. Real physical 
quantities have a subscript r and the constant units have a subscript u. The 
relation between a real quantity q T , the related dimensionless quantity q and 
its scale unit q u is then 

0r = (2) 

To enhance the legibility of the equations subscripts of variables that were 
used in the previous papers have become superscripts. The physical constants 
used in the present paper are listed in Table 1. The mean molecular weight p 
is 0.5 with the assumption that the gas is a fully ionized pure hydrogen gas. 
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The total energy density e r is given by 


1 2 , Pr 

e ' = 2 Pr ' + 


( 3 ) 


where p r is the mass density, v r the velocity and p r the gas pressure. Applying 
Eq. (2) yields 

= 7}PuP v\v 2 + (4) 

The dimensionless equivalent of Eq. (3) has the same form as Eq. (3) itself, 


1 2 P 

c=-pv + 


2 1 7-1 

This gives the relation between the units in Eq. (4) 


( 5 ) 


2 

C u — p u — p u v u . 


( 6 ) 


This simple procedure can now be applied to the other equations. The equation 
of state is 


llT r 


Pr = Pr 

P 

where T r is the temperature. This leads to the dimensionless equation 

p= pT. 

With Eq. (6) the temperature unit must be 

T = —v 2 

u n u ' 

The continuity equation in spherical symmetry is 

d , 1 d ( 2 \ n 

~ Pr + —— ( X rPrV r) = 0 . 


( 7 ) 


( 8 ) 


( 9 ) 


(10) 


dt/ T Xldx r 

where t r is the time and x r the distance measured from the stellar centre. With 
the time unit 


t = — 

4 u — 

the dimensionless continuity equation becomes 


( 11 ) 
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The equation of motion in spherical symmetry is 

d . . 1 d / 2/ QM. 2p r 

w, M + ^w, ( x ’ (PrV ’ + *>) = + ~ 


(13) 


The term £/M» is eliminated with the choice 
Then the dimensionless equation of motion becomes 


(*V J + P)) = 


The energy equation in spherical symmetry is 

d 1 d / , , N \ GM-. xv 

dt r Cr + x? £ Jr ( X,tv(£r + Pr) ) “ ~ prVr ~ + ^ 

The quantity Q* 0 * represents the sum of all external energy source terms 

Qtot _ qcond + Qmtch + qrad + qop 


\tot 

f T * 


(14) 


(15) 


(16) 


(IT) 


with Q c ™ d the thermal conduction, Q™ eeh the mechanical heating, Q T r ad the 
radiation losses and Q° ? the opacity heating. Note that the enthalpy flux is 
not an external energy source term, but implicitly included in the lefthand 
side of Eq. (16). The unit for the energy source terms must be 


tot _ PuK 




The dimensionless form of the energy equation is then 


d 


1 d 




dt ' x*dx 


(18) 


(19) 


So far all the dimensionless equations are independent of the stellar parameters. 
But only if each component of Q tot is independent of the stellar parameters 
will this also hold for the energy equation. 

The thermal conduction is 


Q‘ md i — ( j‘ K T sn —'\ 


( 20 ) 
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Thermal conduction can be defined independently of the stellar parameters 
only by choosing the density unit 


p. = KMKffivt/z,. 


The dimensionless form of thermal conduction then becomes 

>dT\ 


Q 


cond 


--±±( 
* 2 dx y 


,2j’5/2 


dx 


(21) 


( 22 ) 


The mechanical heating is given by 

Qmech 


1 6H r 
x 2 dx T 


(23) 


where H r = F r x 2 is conserved in spherical symmetry if there is no dissipation, 
and F r is the local flux of mechanical energy. We choose the unit for H T in 
such a way that the flux F^ ech deep in the photosphere — an input parameter 
for the calculations — is unitv. Hence 


H u = F? teh Rl 

Thus a dimensionless mechanical flux / is introduced, defined by 

f Hu 
PuVtxl' 

With this definition the dimensionless form of the mechanical heating becomes 


(24) 


(25) 


)mcch 


/ dH 


(26) 


^ x 2 dx 

The weak shock-wave dissipation law for the mechanical heating is given by 

l 

(27) 


dH, ( 12H; y 

— = -' h+1) teJ ' 


dx T 

where v is the frequency of the sawtooth waves. The adiabatic speed of sound 
is defined by 

Pr = Prclh, (28) 

which yields the dimensionless form 


P = pc*h 


(29) 
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if the unit for the adiabatic speed of sound is chosen as 

c u = tv (30) 

The dimensionless form of Eq. (27) becomes 



where a is given by 

o = h + l)(— )’>'-■ (32) 

\ 7 / v u 

Both / and a are input parameters because F^ ech and v are physical param- 
eters. For constant / and a from one star to another an identical solution is 
found, and hence a scaling law is found for F^ eeh and v. 

The radiation losses are 

QT*d = ptpnifo)' (33) 

A good approximation for the radiation loss function P™ d for a hot tenuous 
plasma is a T'^Maw. This approximation will be discussed further in Sect. 3. 
A modified radiation loss function p r r ad (T r ) is defined by 

Pt ad (T r ) = p T T ad (T T ) T -1/2 . (34) 

The function p r T ad is a correction term to the T'^Maw approximation for the 
radiation loss function. If its unit is defined as 

= WlifK;' (35) 

the dimensionless form of the radiation losses Eq. (33) is 

Qtad = p 2 t -\/7 p T*i( Tr )' ( 36 ) 

Q rad is still dependent on stellar parameters through the dependence of p rad 
on T r , but only weakly. 

Finally the opacity heating is 


Q? = W d (x r )plP?(T r -,T tff ). 


(37) 
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The dilution factor W d is a function only of x r / R ., so the equality 

W d (x) = W d (x T ) ( 38 ) 

holds when the distance unit is chosen to be proportional to R+. This is also 
demanded by the boundary conditions at the photosphere, which can only be 
defined independently of the stellar radius, if the distance unit is proportional 
to it 1 . With the definition of the unit for the opacity heating function 

p ? = ( k / h ? A '; 1 t :' /2 ( 39 ) 

the dimensionless form of the opacity heating becomes 

Q 0 * = \Y d {x)p 7 P 0 *(T r -T tB ). ( 40 ) 

The opacity heating P° v depends on stellar parameters through its dependence 
on T r and T tS . In contrast to the radiation losses this dependence is rather 
strong. This is discussed in the next section. 


The complete set of dimensionless equations that describe a one-dimensional 
spherically symmetric stellar atmosphere heated by weak shock-wave dissipa- 
tion is given by 


1 , P 

( 0 ) 

“O 

II 

-O 

(*) 


(c) 

I ( H + ^£(^W+p)) = -^ + f 

(d) 


(e) 


= -hi (*’ rS/ ’f ) ~hf- + w ^ p ° p 



(/) 

(9) 

( 41 ) 


JThis is in contrast to Hammer (1984) who choses the distance unit x u equal to R». This 
is an unnecessary constraint. In the present paper the distance unit x u is chosen to be 10 
Rt rather than J2», see Table 2 
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The units, the relation between the units, and the scaling laws for the physical 
variables are given in Table 2. The unit for the distance is chosen to be 10 R, 
instead of R+ because this gives typical coronal values for the other units. The 
scaling laws are given for two cases, firstly for the general case when the results 
are a good approximation and secondly for a constant A//i?-ratio, in which 
case the results are exact. 


3 Discussion 


The dimensionless equations (41a-d) are independent of the stellar parameters. 
As far as those equations are concerned, the results for one star can be scaled 
exactly to any other star. The energy source equations (41e-g) contain terms 
dependent on the stellar parameters. These are the parameters f and a. 
the modified radiation loss function p rad (T r ) and the opacity heating function 
P^iTr’^Tts). When the free parameters / (dimensionless mechanical flux) and 
q (dimensionless frequency of the sawtooth waves) are constant from one star 
to another by a proper scaling of the related free parameters F™ ech and v. the 
results of a numerical calculation apply to all those stars. 

Expressed in terms of stellar parameters the dimensionless mechanical flux / 
is 

Fr^ ( io^ y /2 (MR^ 

1 10K V hGMq ) (M./M e ) 7/! ' 1 

The factor 10 is the proportionality constant between x u and R* (cf. Table 2). 
Numerically / is given by 

f = 3.91 10-* F? tdl (43) 

which yields / = 100 for the hot 0*star ( Ori used by Korevaar (1989). To keep 
/ constant during the scaling from one star to another, the initial mechanical 
flux F^ eeh must scale like any other energy flux (cf. Table 2) 


F? teh ~ M 7/2 R~ 9/2 . 


(44) 
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The parameter a can be regarded as a dimensionless reciprocal scale length for 
dissipation. A larger value for a gives a larger dissipation and thus a shorter 
effective dissipation scale length (cf. Eq. (41f )). Expressed in terms of stellar 
parameters a is 

<45, 

7 GM e J (A/*/M©) 1/2 

Numerically this is 

o = 3.6 10 s f 11 v (46) 

With v = 1/17000 Hz for £ Ori, a is 6.2 10 3 . To keep a constant during the 
scaling from one star to another, the frequency of the waves v must scale as 

v ~ M 1/2 R~ 3/2 . (47) 

For the rest of this paper it is assumed that F 0 mec/l and v scale as in Eqs. (44) 
and (47). This means that the dimensionless parameters / and a remain 
constant during the scaling, and the results of Korevaar (1989) can be applied 
directly to different objects. 


Given this assumption only the modified radiation loss function p rad (T r ) (Eq. 
(34)) and the opacity heating function P op (T r ;T t ff) (Eq. (37)) depend on the 
stellar parameters since they cannot be defined as function of the dimensionless 
temperature T independently of the choice of the temperature unit T u . 


Opacity heating is completely negligible in a typical corona. It is only im- 
portant for high density — low temperature conditions which are found in the 
photosphere and inside the density pulse that falls to the star during the col- 
lapse phase in the relaxation oscillation. In these cases the gas is in radiative 
equilibrium. The opacity heating then defines the radiative equilibrium tem- 
perature T re in combination with the radiation losses, such that 

P?(T„-,T, t ) = P:“(T„). (48) 

The equilibrium temperature T rt is roughly proportional to the effective tem- 
perature of the star. 
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Consider two stars with the same mass and radius, but with different effec- 
tive temperatures. The photospheres of both stars have different scale heights, 
but the stationary coronae of both stars are exactly the same. The coronal 
structure is completely unaffected by small changes in the photospheric tem- 
perature at the base of the transition region. During the collapse phase in the 
relaxation oscillation all the details are the same except that the temperature 
inside the infalling density pulse is different. This does not influence the relax- 
ation oscillation. Whether the temperature inside the infalling density pulse 
is 5000 K or 15000 K is unimportant. The collapse occurs because the density 
pulse is cool and dense compared with the surrounding gas. The density pulse 
falls almost freely to the star. Consequently the details of the opacity heating 
do not affect the validity of the scaling laws. 

The radiation losses are defined independently of the stellar parameters except 
for the correction function p rad (T T ) (cf. Eq. (34)). A 7’“ 1 ^ J -law for the radiation 
losses has been found to be a good approximation in the range of coronal 
temperatures from 10 s to 10' K (Rosner et al., 1978; Endler et al., 1979). 
Therefore the scaling of the results from one star to another will be a good 
approximation. 

There is a class of stars for which the seeding is exact. When two stars have 
the same A//i?-ratio the temperature unit is the same for both stars as Table 2 
shows. This means that they have the same modified radiation loss function 
p rad (T r ). Then the precise form of the radiation loss function is no longer 
relevant and the scaling is exact. 


4 Applications 

The global coronal relaxation oscillation calculated by Korevaar (1989) is but 
one solution carefully selected to obtain a solution lying just above the thresh- 
old of instability. The stable solution and the solution showing the global 
relaxation oscillation can be scaled using the results of the present paper to 
other stars. No results are available yet of coronae heated by still higher fluxes 
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of mechanical energy. Tables 3 and 4 show the results for a corona round 
£ Ori scaled to an F5-type star, a solar coronal hole, an MO supergiant and 
an elliptic galaxy. Table 3 gives the results for the stable stationary extended 
corona obtained with a mechanical flux a few percent below the limit of sta- 
bility, Table 4 gives the results for the relaxation oscillation obtained with a 
mechanical flux a few percent above the limit of stability. 

Main sequence stars have nearly the same M / R-ratio when M is reduced for 
the outward acceleration due to electron scattering. This can be seen imme- 
diately from Table 3 for the hot O-star ( Ori, an F5-star and the sun. Stars 
earlier than AO have radiatively driven winds, while stars later than FO have so 
much magnetic activity, that acoustic heating becomes relatively unimportant, 
except perhaps in coronal holes. F-stars are therefore the best candidates for 
finding the global relaxation oscillations caused by a high flux of mechanical 
energy. For the F-type star the period of the relaxation oscillation will last 
about 3 months, during most of which the corona will be semi-stationary. The 
collapse phase will last about 1 week, during which time the temperature is 
sufficiently high (about 6 10 6 K) to be easily observable by the X-ray emission. 

The solar corona is heated mainly by dissipation of magnetic energy. For a 
coronal hole it is not clear what the exact heating mechanism is. The dissipa- 
tion scale length for Alfven wave heating is much larger than the observed 0.5 
solar radii (Withbroe, 1988). Acoustic heating may therefore play a role in the 
heating of solar coronal holes. Table 3 shows the main properties. The flux of 
4.0 10 5 erg cm~ 2 s _1 at the transition region is close to the flux of 5 ± 1 10 5 
erg cm _2 s~ 1 deduced from the observations by Withbroe (1988). The period 
of the sawtooth waves of 580 s found from the scaling relation (cf. Table 3) 
is also not inconsistent with the idea of sound wave heating of coronal holes. 
The observed dissipation scale length, which is typically 0.5 R® for a coronal 
hole, (Withbroe, 1988), is roughly 5 times larger than the one deduced here for 
purely acoustic heating, and the observed mass loss is 7 times smaller than in- 
dicated in Table 3. Relaxation oscillations have not been observed in the solar 
corona. Probably the heating of a solar coronal hole is not purely acoustic, or 
perhaps the weak shock-wave dissipation used for the calculations of Korevaar 
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(1989) is not applicable. 

Supergiants usually have a large mass loss. The maximum mass loss calculated 
for the MO-supergiant (cf. Table 3) is extremely low, which confirms that 
in winds from supergiants other mechanisms must be important. Radiation 
pressure on the gas itself does not work for cool giants. Radial oscillations 
and radiation pressure on grains in the outer winds are possible candidates for 
wind mechanisms in cool giants. 

Scaling of stellar results to galaxies is in principle possible, since a galactic 
corona is often treated like a pure gas, i.e. a Maxwellian velocity distribution 
is assumed, so that the equation of state for an ideal gas can be used. The 
geometry is in general non*spherical, but when the scaling is restricted to 
elliptic galaxies the spherically symmetric geometry is a good approximation. 
Table 3 shows that for an elliptic galaxy the M/R - ratio is roughly the same as 
for main sequence stars. The typical coronal temperature of a galactic corona 
is thus the same as for stellar coronae, about 10 6 K, and so is the velocity. 
However, the period of the oscillations is far too large (more than the Hubble 
time) to make them observable. Therefore the coronal thermal instability of 
Hearn et al. (1983) and Korevaar (1989) is not observable in the coronae of 
elliptical galaxies and probably not in other galaxies. 

We conclude that a (re)examination of the X-ray observations of F-stars is 
called for at this time to search for the existence of the coronal relaxation 
oscillation. 
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Table 1: Physical constants 


physical constant 

symbol 

value 

unit 

ratio of specific heats 

7 

5/3 

— 

mean molecular weight 

V 

0.5 

— 

gravitational constant 

Q 

6.67 10- 8 

cm 3 s _2 g _1 

gas constant 

n 

8.31 10 7 

cm 2 s -2 K -1 

conductivity coefficient 

I<0 

1.10 io- 6 

g cm s _3 K“ 7 / 2 

solar mass 

M& 

1.99 10 33 

g 

solar radius 

Re 

6.96 10 10 

cm 
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Table 2: Scaling laws. M and R are the ratio’s between two arbitrary stellar 
masses and radii respectively. 


physical variable 

scaling law 

name 

unit 

relation 

general 

M = R 

distance 

X ^ 

10 

R 

R 

velocity 

t>U 

(GMJx u ) 1 ' 2 

M'f*R ~ 1 ' 2 

1 

speed of sound 

C u 


M i/2 R~ l/2 

1 

temperature 

T u 

WR>) * 

MR - 1 

1 

time 

K 


M-1/2&/2 

R 

density 

Pu 

K 0 (p/n) V2 </*u 

M 2 R~ 3 

R- 1 

pressure 

Pu 

Pu^i 

M 3 R~ 4 

R~ l 

energy density 


Pu 

M 3 R~ A 

R- 1 

energy flux 


Puii 

M 7/2 R- 9/2 

R- 1 

mech. flux x 2 

H . 


M 7 ' 2 R~ 5 ' 2 

R 

radiation loss f. 

P? d 

m/p))VJ<o 

1 

1 

opacity heating f. 

pop 

u 

-.rad'p-1/2 
ru ti 

M-'' 2 R'* 2 

1 

mass loss rate 

M u 

4 TpuVutl 

M*' 2 R~ 3 ' 2 

/r 1 
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Table 3: Scaled properties for maximal stable coronae obtained with a flux of 
mechanical energy just below the limit of stability. The stellar parameters are 
taken from Allen (1973). R. is the radius, M, the mass reduced for electron 
scattering, T max the maximum temperature, M the maximum mass loss, F™ tch 
the mechanical flux at the transition region and P = \/v the period of the 
sawtooth waves that heat the corona by dissipation. 


object 

R, 

(Re) 

K 

(M s ) 

j'max 

(A') 

M max 

(M 0 /yr) 

pmtch 

(erg cm~ 2 s~') 

P 

(s) 

C Ori (0 9.5) 

27.8 

25.1 

2.0 10 6 

1.5 10" 12 

1.0 10 4 

1.7 10 4 

F5-star 

1.2 

1.3 

2.4 10 6 

1.0 10~ 13 

4.4 10 5 

6.7 10 2 

coronal hole 

1 

1 

2.2 10 6 

'T 

7 

o 

o 

4.0 10 5 

5.8 10 2 

MO-supergiant 

12.4 

500 

5.5 10“ 

3.4 10‘ 15 

1.9 10" 3 

1.8 10 6 

Elliptic galaxy 

~10 n 

~10 n 

1.1 10* 

2.5 10~ 3 

1.8 IQ" 7 

1.6 10 14 
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Table 4: Scaled properties for the global relaxation oscillation obtained with a 
flux of mechanical energy just above the limit of stability. t osc is the averaged 
global oscillation period. 7^“ T is the maximum temperature measured dur- 
ing the collapse of the corona. T™° T is the maximum temperature measured 
during the rebuilding of the corona. is the maximum inward velocity mea- 
sured during the collapse and t>™ x is the maximum outward velocity measured 
during the rebuilding of the corona. 


object 

t-oSC 

(yr) 

rpmQT 

1 col 

(A') 

'T’mox 

J reb 

(A') 

C" 

(km/ s) 

,,max 

V Ttb 

(km/ s) 

c Ori (0 9.5) 

7.0 

5.0 10 6 

7.0 10 6 

-500 

700 

F5-star 

0.28 

6.0 10 6 

8.4 10 6 

-550 

770 

solar coronal hole 

0.24 

5.5 10 6 

7.7 10 6 

-525 

735 

MO-supergiant 

760 

1.4 10 5 

1.9 10 5 

-85 

115 

Elliptic galaxy 

2.2 10 12 

2.8 10 6 

3.9 10 6 

-370 

520 
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Abstract 

The current sheet at the top of the postflare loops in a solar two-ribbon flare is shown to bear 
a strong resemblance to the geomagnetic tail. Consequently, in analogy with the geomagnetic 
tail, the generation of high energy neutral particle beams is expected, that leave the sheet with 
a very 6 mall pitch angle. The collisionless direct electric field acceleration of the protons and 
electrons is limited only by the gyromotion of the particles along the magnetic field component 
perpendicular to the sheet. 

Approximately 10 35 protons s ” 1 are produced, and numerical results show a typical energy 
of about 200 KeV, and a power law energy distribution up to 20 GeV. This population is capable 
of generating the observed Hard X-rays and 7 -rays upon impact on the chromosphere. 


Introduction 

Two- ribbon flares are the largest flares observed on the sun, and filament eruptions are often 
considered as their physical cause (Carmichael 1964, Van Tend and Kuperus 1978). Observations 
(Dennis 1985, Cliver el al. 1986), as well as numerical models (Kuin and Martens 1986, Mikic 
el al. 1988), and theoretical considerations (Pneumann 1980, Kaastra 1985) demonstrate the 
origin of a current sheet below the erupting filament, where the postflare loops are formed, and 
where most of the conversion of magnetic energy to particle kinetic energy takes place. This 
current sheet will be considered here in detail. 

During a large two-ribbon flare about 10 32 erg of magnetic energy is dissipated and eventually 
radiated away over the whole electromagnetic spectrum, with the bulk of the emission in X- 
rays. The peak conversion rate, in the impulsive phase of the flare, is about 10 29 erg s -1 . 
The kinetic energy of the ejected filament is of the 6 ame order as the emission, underscoring 
the physical connection between the flare and the filament eruption. A synopsis of the major 
physical processes during a two-ribbon flare is given in Fig. 1 (from Martens and Kuin, 1988). 
The electric, field in the sheet is perpendicular to the plane of the figure below the filament. The 
width of the sheet is in the vertical direction, and the sheet's thickness (negligible compared with 
the dimensions of the flare) is in the horizontal direction. 

The Particle Acceleration Mechanism 

To explain the observed Hard X-ray emission from flares as the result of the impact of a particle 
beam on the dense chromosphere, the mean energy of the particles has to be about 30 KeV if the 


beam is composed of electrons (Lin and Schwartz 1987), and 200 KeV if protons are the carriers 
of the beam energy (Simnett 1986). After acceleration in the current sheet, the mean free path 
(m.f.p.) of electrons as well as protons would be about 10 n cm, larger than the whole erupting 
structure. (Explicit expressions for the calculation of the m.f.p., number of accelerated particles, 
and other quantities discussed below are given in Martens 1988, and shall not be repeated here.) 
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Figure 1. The structure of a two-ribbon flare and the location of the major energy conversion 
processes, viewed in a cross section perpendicular to the filament axis . 


Given the particle energies, and an estimate of the electric field strength in the sheet of 2 
Volt/cm, found from the observed photospheric field strength and ribbon velocities with the use 
of flux conservation (Kopp and Poletto 1986), one can easily calculate the acceleration length 
of the particles in the sheet, by equating the required particle energy to the potential drop in 
the sheet over an acceleration length. The result for the solar sheet is an acceleration length of 
the order of 10 5 cm, 6 orders of magnitude smaller than the m.f.p.! The actual m.f.p. of the 



particles may be 2 or 3 orders of magnitude smaller than 10 n cm, because of collective (e.g. 
wave-particle) interactions, giving rise to anomalous resistivity, but an increase of 6 orders is 
unlikely for anomalous processes. 

* 

Clearly then, particles do not interact during their acceleration process, and one has to 
consider a different way to limit the acceleration path to about 10 s cm. This situation is quite 
similar to that in the geomagnetic tail, where also the m.f.p.’s for protons and electrons are many 
orders of magnitude larger than the tail’s dimensions (see Lyons and Williams 1984 for a review). 
Early model calculations (Speiser 1965, 1967) prompted by in-situ satellite measurements of 
electric and magnetic fields have shown that in the geotail the orbits of individual particles are 
determined by the component of the magnetic field perpendicular to the sheet (B p ): after about 
6ne gyration about this component the particles are ejected from the sheet with a very small 
pitch angle. Protons and electrons acquire the same velocity, and the beams emanating from 
the sheet are consequently neutral. This neutrality, and the small pitch angle, alleviate most of 
the problems connected to the propagation of pure electron beams (see Dennis 1988 for a recent 
review), making proton beams an acceptable candidate. Since the protons and electrons in the 
beam have the same velocity, the protons carry the bulk of the kinetic energy, and the term 
proton beams is appropriate. 


It was shown by Martens (1987, 1988) that the results of the numerical calculations of Speiser 
(1965), and in particular the approximate analytical expressions derived there, can be applied 
to the solar case. The results of this analysis are: 

1. The convection of plasma into the sheet continuously supplies the needed 4xl0 35 protons/sec. 

2. The sheet magnetic field (about 200 Gauss) is acceptable despite the large number of accel- 
erated protons, in contrast with Spicer’s (1983) conclusion. 

3. Neutral particle beams are injected at the top of the postflare loops with almost zero pitch 
angle. 

4. The generation, timing, and spatial distribution of 7-rays and Hard X-rays, can be explained 
from a single power law proton spectrum, with cut-offs at 100 KeV and 20 GeV. 

5. The effective resistivity in the current sheet is over 6 orders of magnitude larger than the 
classical Spitzer resistivity, making the total sheet resistance derived required by observations 
consistent with the microscopic physics. 

Recently we have derived some new results that strengthen the conclusions above. The 
thickness of the current sheet (6) may be estimated from the induction law with the assumption 
of a quasistationary sheet, which is based on the observation that the duration of a flare is much 
longer than the Alfven travel time through the sheet. The result is 



430cm. 


( 1 ) 


With the effective resistivity and the plasma drift velocity from Table 1 we find the result of 
Eq. ( 1). The net energy conversion in the sheet (P) can be derived in the same manner as in 
Kaastra (1985). The result is 

P — IEL « 6xl0 28 ery/sec, (2) 


where I is the current in the sheet, E the electric field, and L the length of the sheet. This result 


is in excellent agreement with the observed Hard X-ray emission in the impulsive phase in large 
flares. Note that the result of Eq. ( 2) is independent of the effective resistivity or the sheet 
thickness. Finally, from magnetic flux conservation we derive an estimate for the perpendicular 
magnetic field component. Equating the inflow of magnetic flux to the outflow, one finds 

VinB ghett = Vout Bp . (3) 

The inflow velocity is equal to the drift velocity, but the outflow velocity is not known a priori. 
Its maximum value is the local Alfven speed. For the geomagnetic tail all the quantities in Eq. 
( 3), except V^t have been measured (Lyons and Williams 1984), and V ou t is found to be equal 
to the Alfven speed. This corresponds to maximum energy conversion in the sheet. We assume 
that the same situation holds for the flare current sheet, and the result is B p « 0.7 Gauss. Given 
this value the average energy of the accelerated protons can be calculated with the approximate 
formula given by Speiser (1965), 

U 9 « 2 m p c 2 (E/B p ) 2 » 174 KeV. (4) 

This result is in excellent agreement with the average energy of 200 KeV inferred from obser- 
vations (Simnett 1986), indicating that the assumption V 4 = V^t is correct. Note that the B p 
calculated is the one at the edge of the sheet. Near the center B p is much smaller, even reversing 
sign at the axis. Consequently the protons accelerated near the center of the sheet can acquire 
much higher energies, theoretically up to the full potential drop along the sheet of 20 GeV. 


Numerical Simulations 


The current sheets of the geotail and a two-ribbon flare are similar, but not identical. So far 
we have used the results of the numerical calculations for the geomagnetic tail done by Speiser 
(1965), and extrapolated them for the solar current sheet. The validity of this approach was 
checked by numerically calculating the orbits of the particles for the flare current sheet. The 
parameters used in these calculations are given in Table 1. The numerical code evaluates the 
orbit of a particle in a stationary electric and magnetic field. The equation of motion for the 


protons is 




df 


= e(E + 


VxS 


)• 


( 5 ) 


dt c 

The integration procedure was tested by a comparison with Speiser ’s (1965) results (see Table 
1 for the parameters) and the correspondence was found to be excellent. For the flare sheet we 
find results that are indeed qualitatively very similar to those for the geotail. A typical proton 
orbit is shown in Fig: 2. The particle initially gyrates along a field line. As it reaches the 
6heet it is accelerated and at the same time is slowly deflected by the perpendicular component 
of the sheet magnetic field. The particle leaves the sheet in a direction almost parallel to the 
magnetic field. The approximation formula for the particle energy (Eq. ( 4)) was found to be 
correct within 30%. It was also found that the initial energy distribution of the protons is of 
little importance for the final results, in accordance with Speiser’s results. In future calculations 
we intend to take into account the presence of a magnetic field component parallel to the sheet, 
and the variation of B p along the width of the sheet. 



Table 1: A Comparison between the Plasma Physics of the 
Geomagnetic Tail and the Current Sheet in a Two- Ribbon Flare 


Physical Quantity 


Geomagnetic Tail 


Observations 


2x10 4 Gauss (1,2) 


2.5x10 Volt /cm (1) 


2-100 KeV (1,6) 
(power law) 


2.6xl0 lu cm (8) 


1.3 — 13xl0 10 cm (8) 


0.02 - 5 cm" 3 (11) 


1.5x10° K (13) 


Two-Ribbon Flare 


Sheet Magnetic field 


Sheet Electrical Field 


Energy of Accelerated 
Protons 


Sheet Length 


Sheet Width 
(^ thickness!) 


Density of Ambient 
Plasma 


Temperature of 
Ambient Plasma 


200 Gauss (3) 

(typical for active region) 


2 Volt/cm (4,5) 


200 Kev - MeV range (7) 
(inferred from Hard X-Rays) 


1.2xl0 10 cm (9) 

(inferred from filament length) 


10 s cm (9,10) 

(inferred from width ribbons) 


10™ cm" 3 (12) 

(typical for active region) 


2.0 — 3.0x10° K (12) 
(typical for active region) 


Plasma Drift Velocity 


Mean Free Path Protons: 
before Acceleration 
after Acceleration 


Debye Length 


Acceleration Length 


Effective Resistivity 


Sheet Thickness (Eq. (1)) 


Number of Accelerated 
Protons 


Net Current 


Net Energy Conversion 
(Eq.(2))' 


Alfven Velocity 


Perpendicular Magnetic 
Field Component (Eq. ( 3)) 


Mean Proton Energy after 
Acceleration (Eq. ( 4)) 


Dreicer Electrical Field 


Spitzer Colliskwiai 
Resistivity 


Larmor Radius upon 
entering: Protons 
Electrons 


Proton Velocity 
after Acceleration 


Gyrofrequency: Protons 
Electrons 


Plasma Frequency 


Calculated 


12.5 km/sec 


3.2xl0 15 cm 
7.6xl0 17 cm 



15 km/sec 


5.2x10 s cm 
3.0xl0 n cm 


9.8xl0 -2 cm 


1.0x10° cm 


4.5xl0“ 12 sei 


430 cm 


3.7xl0 16 Statamp | 1.4xl0 21 Statamp 

(consistent with field strength around the sheet) 


8.1xl0 18 erg/sec | l.lxlO 29 erg/sec 

(excellent agreement with observations) 


4.9xl0 7 cm/sec I 4.4x10“ cm/sec 


1.0x10“° Gauss 0.69 Gauss 


1.3 KeV | 174 KeV 

(excellent agreement with observations) 


2.1xlG“ 14 Volt/cxn T 3.0xl0~ 4 Volt/cm 


1.2xl0- ls sec 1.5xl0 -1 “ sec 


5.8xl0 6 cm 
1.4xl0 5 cm 


4.8xl0 7 cm/sec 


0.31 Hz 
560 Hz 


8.2 cm 
0.19 cm 


6.2x10® cm /sec 


3.0x10° Hz 
5.6x10® Hz 


9.0x10“ Hz 



























































Table 1 (continued) 


Physical Quantity 

Geomagnetic Tail 

Two-Ribbon Flare 

Dimensionless Quantities 

Number of Particles 
in Debye Sphere 

4.0xl0 12 

3.9*10 7 

Plasma Beta (outside sheet) 

5.2*10" 3 

4.3xl0" 3 

Plasma Freq. /Gyrofreq. 

9 

1.6 

Sheet Thickness/Proton 
Gyroradius before Acceleration 

34 

52 

Sheet Thickness/Proton 
Gyroradius after Acceleration 

0.88 

1.9 

Electric Field/Dreicer Field 

1.2*10® 

6.7*10 3 

Parallel/Perpendicular 
Compononent Magnetic Field 

20 

290 

Effective Resistivity/ 
Spitzer Resistivity 


3.0*10 6 

Observed/ Calculated Proton 
Energy after Acceleration 


Hi 
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Figure 2. A Typical Proton Orbit in the Current Sheet of a Two-Ribbon Flare. Inside the 
sheet the perpendicular coordinate has been scaled up by a factor 100. 


Conclusions 

The comparison of Table 1 shows dearly that the plasma physical parameters of the sheets 
considered in this paper differ widely. However, the dimensionless quantities calculated in the 
last part of Table 1 demonstrate that the physical character of the current sheets may be very 
similar. In particular the electric field in both sheets is many orders of magnitude larger than the 
Dreicer field, in sharp contrast the situation in the much denser laboratory plasmas (cf. Stenzel 
et al. 1985). Further, the effective resistivity caused by the ejection of highly energetic particles 
from the sheet, due to the presence of a perpendicular component of the magnetic field, is many 
orders of magnitude larger than the dassical Spitzer resistivity, and also much larger than the 
resistivities that can be obtained through anomalous processes. Note that such a perpendicular 
component of the field is necessary in any reconnecting sheet (cf. Petscheck 1964). 

We condude that a direct observational test to discriminate between electron and proton 
beams in a flare is highly desirable. So far such a direct test has not been proposed. 

Studies of coronal plasma processes at the Harvard-Smithsonian Center for Astrophysics are 
supported by NASA grant NAGW-112. 
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THE MAGNETIC FIELDS ON COOL STARS AND THEIR CORRELATION 
WITH CHROMOSPHERIC AND CORONAL EMISSION 
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ABSTRACT. I discuss the results of recent measurements of magnetic fields on cool 
stars and how these measurements relate to stellar “activity”. Special emphasis is given 
to the correlations between the magnetic field strength, surface filling factor, and the non* 
thermal emission from the hot outer atmospheres of these stars. 

1 INTRODUCTION 

One of the major lessons we have learned from solar research is the critical role of magnetic 
fields in governing the structure and energy balance of the solar atmosphere. Evidence for 
the importance of magnetic fields is widespread. A wide variety of phenomena on the Sun, 
including spots, plages, granules, flares, and loops, are closely associated with magnetic 
fields. These inhomogeneities grow and decay in a repeating 11 year cycle, suggesting the 
continual regeneration of magnetic fields through the operation of a magnetic dynamo. 
And, directly pertinent to this conference, magnetic fields are crucial for the existence 
and stability of the million degree solar corona through their role in the heating and 
confinement of the hot plasma. The dominant heating mechanism for the hot plasmas of 
the chromosphere and transition region (TR) is also likely magnetic in nature. 

By analogy with the Sun, magnetic fields should also be important for other cool 
stars. Indeed, the fraction of the solar surface covered by magnetic structures is relatively 
small (at most, 2 % of the Sun is covered by magnetic fields); many stars show evidence for 
considerably more extensive magnetic “activity* . The fraction of the surfaces of RS CVn 
variables covered by spots, can frequently be 100 times larger than the corresponding filling 
factor of solar umbrae (Vogt 1983). The 11 year cycle of magnetic structures on the Sun 
is also observed in stars: at present some 60% of cool dwarfs have well determined cycles 
based on the modulation of chromospheric Ca n emission (Baliunas 1986). The levels of 
emission from hot plasma in stellar chromospheres and coronae, however, is often orders of 
magnitude larger than the Sun (e.g., Linsky 1985). A thorough study of stellar magnetic 
properties would seem to be in order. 

Unfortunately, the measurement of magnetic fields on solar-like stars is not trivial. 


The most direct means of detecting magnetic fields is through observation of their effect 
on spectral lines. The presence of a magnetic field will, in general, lift the degeneracy 
of certain quantum levels in an atom, separating them slightly in energy (the Zeeman 
effect). The exact number and spacing of the “new” levels (and corresponding allowed 
transitions) depend on the detailed quantum mechanics of the configuration, but for simple 
triplet splitting, the shift in wavelength of the two “new” permitted lines from old one 
is given by AAb(A) = 4.667xlO~ ls g e //A J B, where B is the magnetic field strength (in 
gauss) and g t jj is the effective Lande value of the transition. For typical fields of a few 
kilogauss or less, the Zeeman shift is smaller than the Doppler width of the line and the 
Zeeman transitions overlap, making them difficult to detect. Clever use of the contrasting 
polarization properties of the Zeeman components (the Am = ±1 transitions, called a 
components, are circularly polarized, and the Am = 0 transitions, or x components, are 
linearly polarized), unfortunately, is thwarted by the bipolar nature of solar and stellar 
active regions. If mixed polarities exist within the observed region, a proportionate part 
of the magnetic signal is cancelled, lowering the inferred magnetic flux. In solar work, 
this is not always an important effect, since the distance separating regions of opposite 
polarity is often spatially resolvable. For the completely unresolved surfaces of stars, on 
the other hand, the effect of mixed polarities is ruinous. Almost complete cancellation of 
the circular polarization signal results when the complex, mixed polarities typical on the 
Sun were observed as a point source. It is not suprising then, that very few cool stars show 
any net circular polarization (e.g., Borra, Landstreet, and Mayor 1984). 

2 TECHNIQUES AND PRELIMINARY RESULTS 

So how dots the intrepid astronomer directly measure stellar magnetic fields? Linear 
polarization does not cancel completely in integrated light, but is a second order effect and 
difficult to interpret (Landi Degl ’Innocent i 1982). We are left with the unpolarized line 
profile, with its complex blending of not only a and r components from magnetic regions, 
but contributions from non-magnetic portions of the stellar surface as well. The analysis 
of these various contributions is a difficult problem. Robinson (1980) devised a scheme to 
deconvolve magnetically sensitive (large Lande g) lines from low g lines. The ratio of the 
Fourier transforms of the line profiles is modeled; the period of the ratio is proportional to 
B and the amplitude of the ratio indicates the fraction (f) of the stellar surface occupied by 
magnetic regions. Several improvements on this basic technique have since been developed 
(Marcy 1982; Gray 1984; Saar 1988; Basri and Marcy 1988). 

Preliminary results of stellar magnetic measurements are beginning to emerge. Marcy 
(1984) found that the magnetic parameters change as one moves down the main sequence. 
Gray (1984) found similar magnetic parameters, but did not detect fields on any F dwarfs. 
The most rapidly rotating star in his sample, x 1 UMa, showed substantially larger surface 
averaged field (i.e., the product fB) than the others. Marcy and Bruning (1984) could 
not detect fields on active giants, a fact they ascribe to lower field strengths in these 
stars with lower surface gravities (see also Gondoin, Giampapa, and Bookbinder 1985). 
Consistent with this idea, Saar and Linsky (1986) end Saar, Linsky, and Giampapa (1987) 
reported that among dwarf stars, B increased towards later spectral types and higher 
surface gravities. They also found an increase in magnetic flux with rotation rate (see also 


Linsky and Saar 1987) as first suggested by Gray (1985). Finally, f appears to decrease with 
stellar age (Linsky and Saar 1987), implying that changes in the filling factor dominate the 
time evolution of activity. 

8 MAGNETIC FIELDS AND UPPER ATMOSPHERIC EMISSION 

The apparent connection between magnetic fields and rotation is not surprising, for many 
authors have demonstrated evidence for a connection between rotation and (presumably) 
magnetic-related emission from hot plasmas in stellar chromospheres and coronae (e.g., 
Pallavicini et al. 1981; Mangeney and Praderie 1983; Vilhu 1984; Simon and Fekel 1987). 
The age dependence of chromospheric activity has also been extensively studied (Sku- 
manich 1972; Soderblom 1983; Hartmann and Noyes 1987), and thus the observed decrease 

* 1 neters with age is not unexpected. It is important, however, to deter- 

ctions between the magnetic fields and the outer atmospheric activity 
d the underlying physics of the heating. 
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ning the methods of Saar 1988) with X-ray and Ca II H+K fluxes, 
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lual active regions on the Sun show a similar relation between X-ray 
Schrijver 1987b), suggesting that the coronal heating mechanism is 
letic flux density. Note that since the magnetic analyses methods 
,h considerable magnetic flux, the Sun is an important calibration 
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Figure 1: Soft X-ray flux density F* vs. the mean magnetic flux density fB. The solid 
line is F x oc(fB) 0#9S . Magnetic and X-ray fluxes at solar activity minimum and maximum 
are connected with a dotted line. For comparison, the mean value of solar active regions 
(symbol at fB = 100 G) is shown (Schrijver 1987b). Symbols: circles, squares and triangles 
mark G-, K- and M-type dwarfs respectively. 



Figure 2: Ca II H+K excess flux density AFc<»/j vs. fB. Symbols as in Fig. 1. The left 
segment of the solid line is AT call « (fB) 0 * 62 , while the right segment is an approximation 
to the F* - AT call saturation seen in active stars (Schrijver 1987a). The horizontal line 
segment shows the maximum value observed in solar active regions (Schrijver and Cot£ 
1987) and G dwarfs. The color-dependent decrease in chromospheric emission from M 
dwarfs has been corrected for by the deficiency factor 7 (see Schrijver and Rutten 1987). 




I have also studied correlations between the preliminary magnetic parameters and 
other diagnostics of chromospheric and TR plasma. Stellar surface fluxes in O I (1304 A), 
C II (1335 A), C IV (1550 A), and Si II (1808 A) measured with the IUE satellite and 
analogous solar measurements were gathered from Oranje (1983) and a number of other 
sources. The preliminary results (based on a small sample of stars) are: F civ oc(fB) 0,55 ; 
Ten oc(fB) 0 - 44 ; F 5,7/ oe(fB) 0 - 40 ; F oi oc(fB) 0 - 36 (see Figure 3). The cooler chromospheric 
lines (e.g., O I, with T/ 0 rm » 6x10 s °K) show a reduced dependence on fB relative to lines 
formed higher in the atmosphere (C II, T/ OTm 2xl0 4 °K, and C IV, T f„m « 10 5 °K), 
consistent with the F, and F call results. Note, however, that no “basal” correction has 
been applied to the O I and Si II data. 

Once a larger set of magnetic measurements has been made, it should eventually 
be possible to begin to test various atmospheric MHD wave heating models for chro- 
mospheres, TR and coronae (e.g., Ulmschneider and Stein 1982), test empirical activity 
theories (Vilhu 1984) and predictions of the range of dynamo time behavior (e.g., Noyes, 
Weiss and Vaughan 1984) and guide the development of more detailed theories. 

4 ROTATIONAL MODULATION OF MAGNETIC AND RADIATIVE FLUXES 

Observations of magnetic fields and a wide range of other activity-related diagnostics (Ca 
II, IUE line emission, He I D3, broaadband linear polarization) on the active G8 dwarf, £ 
Boo A, were made in June of 1986. Once again, as no compensation for weak blends was 
made, the resulting magnetic fluxes should be regarded as preliminary. 

Figure 4 shows the normalized magnetic (« fB) and ultraviolet fluxes for £ Boo A 
plotted versus rotational phase (P = 6.2 days, Noyes et al. 1984; <f> = 0 defined as the phase 
of maximum Ca II emission). Rotational modulation of the emission line and magnetic 
fluxes, in phase with one another, is clearly visible. Lines formed at higher temperatures 
(C IV, Si IV, He n, C n) appear to show greater modulation than the cooler lines formed 
in the chromosphere (Ca n, Mg n, Si II, C I, O I). The greater sensitivity of hotter lines 
to magnetic flux in rotationally modulated plages is consistent with the results in §3 for 
entire stars, suggesting that perhaps the main difference between active and quiet stars 
is the number or size of plages, and not a major difference in their physical properties. 
The reduced sensitivity of the chromospheric lines to changes in magnetic flux suggests 
once more that for stars as active as £ Boo A (fB m 600 Gauss), the heating mechanism 
for the chromosphere “saturates” , while the heating mechanisms operating higher in the 
atmosphere still are strongly affected by the variable magnetic flux. 

Further information on the nature of the magnetic regions on £ Boo A can be derived 
from the contemporaneous measurements of broadband linear polarization (Pl)> While 
the polarization amplitude varied by a factor of 3, the angle which the linear polarization 
vector made with respect to the N-S axis never varied from the 0-180 degree line by 
more than about 10 degrees. The simplest configuration that can reproduce this result 
(for fields normal to the surface) is a series of active regions which cross the stellar disk 
near disk center. Since sin t « 0.6 for £ Loo A, the active regions would then be located 
at intermediate latitudes on the average. The line perpendicular to the polarization angle 
(i.e., at « 90°) then corresponds to the stellar rotation axis, entirely consistent with the 



Figure 3: IUE emission line surface flux densities vs. magnetic flux density. The solid lines 
are the power law fits F C /v oc(fB) 0 - 55 , F C // oc(fB) 0 - 44 , and F 01 oc(fB) 0 - 36 . Data for the 
quiet and active Sun are connected by dashed lines. Symbols as in Figure 1. 








Figure 4: The magnetic flux (a fB, heavy solid), various ultraviolet emission line fluxes (C 
IV, solid; C II, dashed; O I, long dashed) and the smoothed, Ca II excess flux (dash-dotted) 
of £ Bootis A, normalized and plotted versus rotational phase. 

known orbital axis (at 78°) of the £ Boo A+B system. 

The amplitude of P i measures the net tangential component of the magnetic field and 
is maximized when significant magnetic flux is distributed asymmetrically, near the stellar 
limb (e.g., Landi Degl’Innocenti 1982). The combination of the variation of fB (weighted 
towards disk- center regions because of projection and limb-darkening effects) with the am- 
plitude of P& therefore contains positional information, which Saar et al. (1987) used to 
construct a crude "magnetic image” of £ Boo A. Using more sophisticated modeling, we 
hope to map the active region geometry of £ Boo, which, when combined with the radiative 
and magnetic fluxes as a function of phase, will allow us for the first time to derive realistic, 
two-component (quiet + active) models for an outer atmosphere of a solar-like star. 
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Abstract 

In this paper the effects of the parameter u> p e /fi e on the electron cyclotron 
maser instabilities are calculated. The range of this parameter is divided 
into several parts to discuss the polarization sense and degree of the solar 
radio spikes. Five possibilities for the harmonic structure of radio spikes are 
predicted (the frequency ratio is about 1:2 or 2:3). The conclusions of this 
paper are compared with the relevant observations. 


1. Introduction 

As a fine structure of solar activity and a direct signature of energetic par- 
ticles in solar flares the observations of solar radio 6pikes have been stud- 
ied by many scientists (Melrose, 1984; Sharma et a/., 1984; Vlahos et al . , 
1985; Li, 1986, 1987; Huang, 1987, 1988a, b). The theory of the electron 
cyclotron maser instability seems to be a reasonable model for explaining 
the millisecond radio spikes in the decimetric range (Benz, 1986). Many 
observational characteristics of radio spikes are consistent with the predic- 
tions of this theory. For instance, the high brightness temperature of the 
radio spikes is connected with the high growth rates of the electron cyclotron 
maser instability , and because the electron cyclotron maser instability is a 
resonant wave-particle interaction in velocity space, the band-width of the 
radio spikes is very narrow, and the directionality of 6pike emissions is very 
strong. 
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Some aspects need further theoretical consideration, such as the polar- 
ization and the harmonic structure of radio spikes. The purpose of this 
paper is to explain these important properties by the electron cyclotron 
maser instability theory. 

In Section 2, the effects of u >pe/Sl e on the electron cyclotron maser in- 
stability are discussed and calculated. The polarization of the electron cy- 
clotron maser instability is discussed in Section 3. The possibilities for the 
harmonic structure in the radio spikes are predicted in Section 4. Section 5 
is a summary. 


2.The Effects of the Parameter ui^/Qe on the electron cy 
clotron maser instability 


The formula for the growth rates of the electron cyclotron maser instability 
was given by Wu (1985) for the hollow beam distribution of the non-thermal 
electrons 
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Here, u>j and u T are the imaginary and the real parts of the frequency, where 
it is assumed that u, < u> r . and Sl e are the electron plasma and cyclotron 
frequencies, respectively. N is the refractive index, k± and are the vertical 
and horizontal wave-vectors, B is the angle between the propagation vector 
and the magnetic field. J n (6 n ) is the usual Bessel function of order n. £o = 
— »'uoxA*_l, 2, is the plasma dispersion function, tio|| and u 0 x are the bulk 
flow momentum parallel and vertical to the magnetic field, aj| and «x denote 
the characteristic momentum spreads. In addition, 
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It is noted that there are two basic parameters of the ambient magnetized 
plasma: the electron density n e and the intensity of the magnetic field Bo. It 
is easy to show that the growth rates of electron cyclotron maser instability 
depend only on the ratio of w,* and fl e , or the ratio of n X J 2 and Bq. 
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From the formula for the cut-off frequencies in a "cold” plasma (Kong 
et al., 1982), 

(1) 2 

= (- ^ + 0.25) 1 / 2 * + 0.5 ( Fast Extraordinary Modes), ( 5 ) 

( 2 ) 

(Ordinary Modes), (6) 

it is evident, that the lower harmonics will be suppressed gradually with 
increasing u)pe/Sl e . 

The question is how to choose the range of Upe/(l e for the calculation of 
the growth rates of the electron cyclotron maser instability . It is supposed 
that solar spikes are mainly composed of the second harmonics of ft e . And 
the intensity of the ambient magnetic field is estimated as 50-500 Gauss 
corresponding to the decimetric range (300-3000 MHz). The average value 
of n e in a solar active region is about 10 9 cm -3 (Mclean, 1986) so the value 
of ufpe/Sle should be in the range of 0.2-2. 

On the other hand, the value of the plasma (3 in the low corona is about 
10~ 4 (Dulk et al., 1978) and 


a 2n e kT e 


ill m e c2 


( 7 ) 


Here k is the Bolzhman constant, T e is the electron temperature, y is the 
magnetic permeability, and c is the velocity of light. 


With T e a 10 6 K, 4kT e /m e c 2 ~ 10~ 4 , and Upe/fle ~ 1» this value is 
consistent with the estimates in the above paragraph. The calculations for 
the growth rates of the electron cyclotron maser instability at the first and 
second harmonics of the fast extraordinary and ordinary modes are shown 
in Figures 1-4. It is evident from these figures that: 


1. All harmonics of all modes can only be excited in a respective range 

of Ljpg/ilf. With increasing u> pe /Sl e , the growth rates of the electron cy- 

clotron maser instability vary non-monotonically (there is one maximum of 

the growth rates for each curved line). 
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2. The range of w pe /fl e for the fundamentals of the fast extraordinary 
modes is from 0.2 to 0.5 (Fig.l); for the second harmonics of the extraor- 
dinary modes this range is from 0.5 to 1.3 (Fig.2); for the fundamentals of 
the ordinary modes from 0.4 to 0.95 (Fig.3); for the second harmonics of the 
ordinary modes from 0.5 to 1.7 (Fig.4). 

3. The range of Wpe/f l e can be divided into several parts: For the values 
between 0.2 and 0.5, there are only the fundamentals of the fast extraordi- 
nary modes to be excited; for the values between 0.5 and 0.9, the fundamen- 
tals of ordinary modes and the second harmonics of the fast extraordinary 
modes are both excited; for the values between 0.9 and 1.2, the growth rates 
of the second harmonics of the fast extraordinary modes are one order of 
magnitude higher than that of the second harmonics of the ordinary modes; 
for the values between 1.2 and 1.7, there are only the second harmonics of 
the ordinary modes to be excited; beyond that, the third harmonics of the 
fast extraordinary modes and the ordinary modes will be excited gradually. 

These results are very useful for the discussion of the polarization and 
the harmonic structure of solar radio spikes. 


3. Polarization 


The polarization of solar radio spikes is observed by Benz (1982, 1987), but 
it is not discussed in theoretical detail. In fact this property can be explained 
directly by the electron cyclotron maser instability theory. 

According to Mdean (1986) the polarization degree is defined as: 


_ Tbr - Tbl 
Tbr + Tbl 


( 8 ) 


Here, p is the polarization degree (%), Tbr is the brightness temperature 
of the right-hand circular components, Tbl * s the brightness temperature of 
the left-hand circular components. 


On the other hand, the solutions of the ’’cold” plasma dispersion equation 
(Wu, 1985) are: 


A 0 (ifc,u>) = AN 4 + BN* + C = 0, 

(9) 

„ r2 -B ± VW* - 4 AC 
K * = 2 A 

(10) 

A = e sin 2 6 + 17 cos 2 0, 

(11) 


ft 
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B = — eri(l + cos 2 0) — (e 2 — g 2 ) sin 2 0, (12) 

C = ife? - g 2 ). (13) 

Here N+ and iV_ are the refractive index of the ordinary and the extraor- 
dinary waves respectively. It is easy to show that the left-hand circular 
(right-hand circular) waves and the ordinary modes (the fast extraordinary 
modes) correspond to the same refractive index. This result is consistent 
with the ‘Leading Spot Rule’ (Benz, 1987). If the growth rate of electron 
cyclotron maser instability is constant during the rise-time of a single spike, 
then the formula for the polarization degree can be written as 

_ Tbx - Tbo 
P ~ Tbx + Tbo ’ 

pUfixAkT _ pVio&T 

= 1 (14) 

Here, Tbx » T BO are the brightness temperature of the extraordinary and the 
ordinary waves, respectively, A r is the duration of the rise time of a single 
spike, and the resonance absorption is neglected in this formula. Because 
the growth rates of the electron cyclotron maser instability are different for 
different values of Upe/Sl e , on the basis of the calculation in Section 2, the 
range of w pe /$l e can be divided into several parts to analyze the polarization 
degree. 

When the value of u> pe /Q. e is from 0.2 to 0.5, there are only the funda- 
mentals of the fast extraordinary modes with right-hand circular polarization 
and the polarization degree p is very high (up to 100%). But, when u>pe/n e is 
larger than 0.4, the growth rates of the fundamentals of the ordinary modes 
will increase gradually, which will result in a decreasing p. 

When the value of u>pe/n e is from 0.5 to 0.95, there are only the funda- 
mentals of the ordinary modes with the left-hand circular polarization and 
the p- value is up to 100%. 

When the value of Wpe/ft* is from 0.5 to 1.3, both the fast extraordinary 
and the ordinary modes at the second harmonic may be excited, but the 
growth rates of the former are four to five times larger than that of the 
latter. If the value of w^/fle is larger than 1.1, the differences between these 
two modes decrease gradually, and the polarization degree of the right-hand 
circular will vary from 75% to 15%. 


f 
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When the value of u pe /Sl e is from 1.3 to 1.7, the second harmonics of 
the fast extraordinary modes will be suppressed rapidly, and the second 
harmonics of the ordinary modes decay slowly. Hence, the polarization 
sense will be changed into the left-hand circular, and the polarization degree 
will increase from 0 to 100% with the decrease of intensity of the ambient 
magnetic field. 

The high circular polarization in some observations of solar radio spikes 
(Benz, 1986) can be explained by the fundamental emmisions of the fast 
extraordinary modes, which propagate at a small angle 6 and the resonant 
absorptions are very weak (Huang, 1987). More recent measurements (Benz 
et a/., 1982; Stahli et al., 1986; Nomino et a/., 1986) agree that the po- 
larization of solar radio spikes can vary from 0 to 100% with left-hand or 
right-hand circular senses at different frequencies, which may be correspond 
to different values of u pe /Sl e in the sources of radio spikes. It is interesting 
to find the polarization sense sometimes changing from left circular to right 
circular (Benz et al., 1982). A possible explanation is that when the electron 
beam propagates along the open magnetic field lines, with the decrease of the 
value of Wpe/fte, the fast extraordinary modes will be suppressed gradually 
instead of the ordinary modes (Figures 1-4). 

On the other hand, from the measurements of the polarization sense and 
degree of the radio spikes, the value of u pe /Q e (also the plasma density) in 
the sources of spikes can be estimated directly. 

4. Harmonic Structure 

The harmonic events of solar radio spikes have been observed by Benz (1987), 
Stahli et al. (1986), and Yun Nan Observatory of the Academia Sinica (Xie 
et al., 1987). However, the theory of the harmonic structure is still unsolved. 
The author indicated (1987), that the first and the second harmonics of the 
fast extraordinary modes cannot be observed in the same direction, because 
their 0- values are quite different (Fig.5). From the calculations in Section 
2, the harmonic structures of radio spikes can be excited in several regions 
of 

When the value of u pe /il e is from 0.53 to 0.92, and the 6 - value is about 
70°, the growth rates of the ordinary modes at the first harmonic and the fast 
extraordinary modes at the second harmonic have similar magnitude. So it is 
possible to observe the harmonic structure of the frequency ratio 1:2. Their 
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polarization senses are opposite, and the amplitude of the fundamentals is 
larger than that of the second harmonics. But, when the value u pe /Cl e is 
larger than 0.8, the amplitude of the second harmonics will exceed that of 
the fundamentals (Fig.6). 

When the value of u)pe/Sl e is from 0.91 to 0.97, and the angle 9 is about 
80°, it is possible to receive both the fundamentals and the second harmonics 
of the ordinary modes with the same polarization sense (left-hand circular), 
and the frequency ratio is about 1:2 (Fig. 7). 

When the value of u>pe/fl e is from 1.3 to 1.35, the extraordinary modes 
at the second and the third harmonics may be observed in the same direc- 
tion (9 is about 80°) with the frequency ratio 2:3 and right-hand circular 
polarization (Fig. 8). 

When the value of is between 1.7 and 1.9, the second and third 

harmonics of the ordinary modes also can be received at an angle 9 of about 
70° (2:3 and left-hand) (Fig. 9). 

When the value of Upe/fte is between 1.5 and 2.0, the ordinary modes 
at the second harmonic and the extraordinary modes at the third harmonic 
may be observed with the opposite polarization sense and the directional 
angle 9 is about 70° (Fig.10). 

There are three reasons why it is not very easy to observe the harmonic 
structure of the solar radio spikes: 1. The range of Wpe/Sl e is very narrow 
for the harmonic events. 2. The range of the directional angle 9 is also very 
narrow (often much smaller than ten degrees). 3. Because the harmonic 
events can only be received at a large value of the angle 9 , the resonance 
absorptions are very strong. 

The properties of the harmonic structures of radio spikes observed by the 
Yun Nan Observatory can be explained by the emissions of the fundamentals 
of the ordinary modes and the second harmonics of the fast extraordinary 
modes, whose frequency ratio, amplitude ratio and directional angle are 
consistent with the observations (Xie, et al. 1987). The harmonic structures 
with the frequency ratio of 1 : 1.39 and left circular polarization (Benz, 1987) 
can be explained by the emmisions of the second and third harmonics of the 
ordinary modes (Figure 9). 
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5. Summary 

In this paper, the effects of the parameter uy/fte on the electron cyclotron 
maser instability are emphasized and calculated, which is the basis of the 
discussion of the polarization and the harmonic structure of the solar radio 
spikes. The polarization properties (including the sense and the degree) are 
analyzed in several regions of fi e . A high right-hand circular polariza- 
tion (up to 100%) exists in the range from 0.2 to 0.5 (the fundamendals of 
the fast extraordinary modes). A high left-hand circular polarization (up to 
100%) exists in the range from 0.5 to 0.95 (the fundamentals of the ordinary 
modes). A middle right-hand circular polarization (75% to 15%) exists in 
the range from 0.5 to 1.3 (the fast extraordinary and the ordinary modes at 
the second harmonics). A middle left-hand circular polarization (0 to 100%) 
exists in the range from 1.3 to 1.7 (the second harmonics of the ordinary 
modes). 

Finally the possibilities for the harmonic structure are predicted in five 
separate regions of Upe/Cl e ' The harmonic events of the frequency ratio 
1:2 may be observed between 0.53 to 0.92 at 8 » 70° with the opposite 
polarization (the left-hand circular for the fundamentals, and the right-hand 
circular for the second harmonics), and between 0.91 and 0.97 with the same 
left- hand circular polarization at 6 « 80°. There are three possibilities for 
the frequency ratio 2:3, the first is 1.3 to 1.35 at 8 « 70° with the same right- 
hand circular polarization; the second is from 1.7 to 1.9 at 9 « 70° with the 
same left-hand circular polarization; the third is from 1.5 to 2.0 at 8 as 70° 
with the opposite polarization sense (the left-hand circular polarization for 
the second harmonics, the right-hand circular polarization is for the third 
harmonics). But, because the ranges of u>pe/ft e and 8 are both limited and 
the resonance absorptions are very strong, the harmonic structures are not 
often observed. 
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Figure Captions 

Figure 1: The relation of growth rates of the fundamentals of the fast 
extraordinary modes and w pe /fi e , 0. («o/c = 0.5, a = 40°, oj. = aj| = 0.2, 
a is a pitch angle between the beam and magnetic fieldln the following 
figures these beam parameters are the same) 

Figure 2: The relation of growth rates of the second harmonics of the 
fast extraordinary modes and u pe /Sl e , 0. 

Figure 3: The relation of growth rates of the fundamentals of the ordi- 
nary modes and Upe/ile , 0. 

Figure 4: The relation of growth rates of the second harmonics of the 
ordinary modes and Upt/Q e , 6. 

Figure 5: The comparison of growth rates of the fundamentals and sec- 
ond harmonics of the fast extraordinary modes. 

Figure 6: The comparison of growth rates of the fundamentals of the 
ordinary modes and the second harmonics of the fast extraordinary modes 
(0 ~ 70°). 

Figure 7: The comparison of growth rates of the fundamentals and sec- 
ond harmonics of the ordinary modes ( 0 « 80°). 

Figure 8: The comparison of growth rates of the second harmonics and 
third harmonics of the fast extraordinary modes ( $ « 70°). 

Figure 9: The comparison of growth rates of the second harmonics and 
third harmonics of the ordinary modes (0 « 70°). 

Figure 10: The comparison of growth rates of the second harmonics of 
the ordinary modes and the third harmonics of the fast extraordinary modes 
(0 ss 70°). 
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Abstract 

The classical theory of the tearing mode instability is discussed and a formula 
for the growth rate of Alfven waves is derived from the dispersion relation for 
a magnetized plasma, in which the plasma resistivity and the electric current 
along the ambient magnetic field are both considered. The growth rates are 
directly proportional to the plasma resistivity, the electric current density 
and a parameter u> pe /ft e . Moreover, the growth of Alfven waves can only 
be excited in a direction nearly parallel to the ambient magnetic field, and 
the velocity of the electron fluid (related to the electric current) should be 
larger than the Alfven velocity. Therefore Alfven waves and tearing modes 
are similar physical phenomena. 


1. Introduction 

The tearing mode instability is one of the most important large-scale MHD 
instabilities. As the prime phenomenon and exciting mechanism in the dis- 
ruptive instabilities in tokamaks, the tearing modes have attracted wide at- 
tention (Kadomtsev, 1984). In plasma astrophysics, the tearing mode may 
also play an important role in the magnetic energy conversion and particle 
acceleration in solar flares (Priest, 1983). 

The classical theory of the tearing mode instabilities was established by 
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Furth et al. in 1963. Because of the shear in the magnetic field and the 
periodic boundary condition, the tearing modes can be excited only in a 
narrow boundary layer of the corresponding rational magnetic surface. So 
the WKB approximation is used: the resistivity of plasma is considered only 
in the boundary layer and the ideal MHD equations remain valid in the other 
regions. Hence the growth rates of the tearing modes are derived from the 
continuity boundary condition. 

The classical theory has been developed further over the years: from the 
plane slab to the cylindrical and toroidal coordinates, from the linear theory 
to the non-linear theory, from the analytical model to numerical calculations. 
However, some questions remain open: 

1) Bateman (1982) derived the same result for the growth rates of the 
tearing modes as in the classical theory (Furth et al., 1963) using a very 
simple model without solving any boundary problem. 

2) In the classical theory, the imaginary parts of the eigenvalues are 
neglected. However, the imaginary parts of eigenvalues describe a periodic 
oscillation of the tearing modes, and the real parts of the eigenvalues describe 
the exponential growth of the amplitude of the oscillation. The periodicity 
of the tearing modes is proven by the author (1987, 1988) from the imag- 
inary parts of the eigenvalues in MHD equations and the real parts of the 
eigenvalues are consistent with the exponential growth rates in the classical 
theory of tearing modes (Furth et al., 1963). 

The quasi-periodicity of the tearing modes was shown directly in the 
disruptions in tokamak plasma (Shunji Tsuji, 1983). The evolution of solar 
radio spikes is also explained by the quasi-periodicity of the induced electric 
field of the tearing modes (Huang, 1988). 

3) The source of free energy for the growth of tearing modes in a toroidal 
plasma is the toroidal electric current along the ambient magnetic field. 
However, this relation is not shown in the formula for the growth rates of 
the tearing modes in the classical theory. 

4) In fact, the resistive boundary layer is not a clear physical boundary 
and its width depends on the properties of tearing modes, which makes the 
growth rates of tearing modes vary with the evolution of the tearing modes 
themselves, so the growth rates in the classical theory of tearing modes 
are not the strict eigenvalues (a strict eigenvalue is a constant without any 
dependence on the solutions). 
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5) Tlie anomalous growth of the tearing modes has not been explained 
very satisfactorily and it is often considered to be the result of the interaction 
between different modes. However, the author has shown that the anomalous 
growth of tearing modes is an inevitable result of the growth of normal 
tearing modes, even for a single mode (Huang, 1987). 

In this paper the author is concerned in the questions 1-4. Question 5 
will be discussed in a forthcoming paper. In Section 2 the growth rates of 
the Alfvcn waves in a current-driven instability are derived from the plasma 
dispersion relation. In Section 3 the growth rates of the Alfven waves are 
compared with the classical theory of the tearing mode instability. Section 
4 is a summary. 


2. Growth Rates of the Alfven Waves 
2.1 Physical Analysis 

Because of the periodicity of the tearing modes, there must be a correspond- 
ing electromagnetic wave in the dispersion relation of the plasma. Note that 
tearing modes have a long wavelength or a low frequency and their phase 
velocity is of the same order of magnitude as the Alfven velocity (Huang, 
1987, 1988). Therefore the growth of the tearing modes may be related to 
the excitation of the Alfven waves, which attract the free energy of the elec- 
tric current of a magnetized plasma. Moreover, the growth of the tearing 
modes only takes place in a boundary layer, which is much smaller than the 
characteristic length of non-uniformity of plasma and much larger than the 
Debye length of plasma. So it is possible to use the method of solving the 
dispersion equation of a uniform plasma. 

Note that there is no dissipation mechanism in the "cold” plasma dis- 
persion relation, so its dispersion tensor is a complex symmetric tensor and 
its eigenvalues will be real. When the resistivity and the electric current of 
the plasma are both considered in the MHD equations, the dispersion tensor 
of the plasma will have an antisymmetric part and its eigenvalues will have 
an imaginary part. 
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2.2 Derivation of Growth Rates 


Initially the effect of ion motions is neglected, and the momentum equation 
for the electron fluid is: 

dV V x5. s 

n e m e — — = —en e (E H ) + R. (1) 

at c 

Here V is the velocity of the electron fluid related the background ions, n e 
is the density of electrons and m e is the mass of an electron. The collisions! 
term R is defined as 

R = en e T)o J (J = -en e V). (2) 


Here q o is the resistivity of plasma. 

Let V = Vo + V\, B = B 0 + B\ and E = Eq + E\ and it is easy to write 
down the expansion of Equation (1) 


e 0+ ^LB, = voJ ^ 


(3) . 


= - e( E, + *1 ‘J° + - e^V,. 

dt c c 


(4) 


Equation (3) is simply the generalized Ohm’s law. Working out Eq. (4) 
yields 

C U r 



When Vo = Voi and k = k x x + k z z, the tensor A can be divided into two 
parts, a symmetric tensor and an antisymmetric tensor 


A = Aj + A2 

l -fe||V 0 0 fc x Vo/2 \ / 0 0 k±V 0 /2 \ 

0 -fc,|V 0 0 + 0 0 0 

\ k ± V 0 /2 0 0 ) \ —k'iVo/2 0 0 / 


(7) 
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Take the Fourier transform of Equation (4) (Melrose, 1982): 


(uSij + tT2 e £,-jj.6fc)Vij _ —ie/m e (6ij + Aij)E\j. (8) 

Here, w = u r + tu»o, u>o = e 2 n e »jo/m e , fi e = Ho = Ho6, the Levi-Civita 
symbol ey* is defined as: 



when ( i,j,k ) is an even permutation 
when ( i,j,k ) is an odd permutation 
when two of i,j, k are of the same value 


Let 


+ iQe£ijkbk)Cjl — U r 6ii, 


the dielectric tensor of the plasma can be derived from (7) 


<*>; 


D ij = 6 ij --%C il (6 lj + A lj ), 


w, 


( 9 ) 



4irn e e 2 

m e 


C = 


(A. °\ 

tti/Qfe co 2 n I 

ZTTUr 0 
0 0 l) 


( 10 ) 


The tensor C can also be divided into a symmetric and an antisymmetric 
part, 

C = Ci + C 2 





0 


— tU/rOe 

if 

^7 

0 



+ 


t'2w T wo 


1 0 0 \ 

0 1 0 I . 
0 0 0 ) 


( 11 ) 


Under the approximations u>, <C u> r , <C (w; and w r are the real and 
the imaginary parts of the frequency, Dfj and Dfj are the antisymmetric 
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and the symmetric components of the dielectric tensor, respectively), it can 
be proven (Melrose, 1982) 


u>, _ — Im (ty • D) 
u> r G 


( 12 ) 


G = 



(13) 


Here Ao = det(D), and ¥ • D = Aol, I is the unit tensor. For the ”cold” 
plasma model (Wu, 1985) 


= N 4 sin 2 6 - N 2 (e sin 2 6 + 77 ) + e*?, 


$ yy = -7V 2 (esin 2 0 + 77 cos 2 0) + erj, 

9 ZZ = N 4 cos 2 6 - N 2 e(l + cos 2 6) + e 2 - g 2 , 

= 2N 2 (N 2 - e)sin0cos0 = 

® y , = i2N 2 g sin 6 cos 0 = — ^l X yi 
Vxy = ilg(V - N 2 sin 2 6) = -9 yxy 

a = - e) 2 ^j + (1 - >))cos J «]+ 

1 1 - 1 ) 

2[e 2 (l + o)-j ! + 22?Ilri£]. 

Here, c = 1 - g = 1 1 ^ » N is the refractive index, 

6 is the angle between the magnetic field and the wave-vector. 

From (9), (7) and (11), it is evident, that Ci • A 2 has no effect on the 
antisymmetric part of the dielectric tensor, so 


u, 


D - s~fCj.(I + Ai), 


u: 


= i\ 


l.Mi 

Ul r 


0 

0 


0 

1 

0 


1 _ 


krVn 

2u *r 
0 


(14) 
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2u> r U 0 u le 

U 2 r~W 


( 15 ) 


Then with (14) and (15) the formula for the growth rates can be written as 


u r 


Al(*~ + - 1) - 

U U) r u r 


(16) 


For Alfven waves N = and u> r < fl e , fii. When 6 = 0 the wave- 

vector is parallel to the magnetic field, and formula (14) can be simplified 
to 


1 / fcgVp _ jx 
U T Tn e L>r Sl\ N 2 (Jj r ' 


(17) 


erjoJo <Yl (1 Y ±' | 

m t u r V A ni c 2 1 V 0 ’ 


(18) 


When 6 = it/2, or fcj| = 0 


fx_ _ e 2 n e 7?o tJp e 1 0 

U> r m e U) T ft 2 


(19) . 


It can be proven that Ui is a monotonic function of the angle 6 between 
0 - it /2. 


3. Conclusions 

1) From Equation (17), when rfy = 0, it follows that u >,• = 0, which means 
that the growth rate of the Alfven waves depends on the resistivity of plasma. 
This result is consistent with the theory of tearing mode instabilities. 

2) From Equation (17), when Vo = 0 or Jo = 0> the growth rates will 
be negative, hence an electric current along the magnetic field is a basic 
condition of the growth for the Alfven waves. The tearing mode instabilities 
also depend on the electric current. 

3) From Equation (16) or (18), the growth rate is positive only when 
Alfven waves propagate nearly parallel to the magnetic field. 

4) From Equation (17), <*>,• a which means that the growth of the 
Alfven waves is suppressed gradually with increasing magnetic field or de- 
creasing plasma density. 
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5) From Equation (17), the Alfven waves of the lower frequencies are 
most easily excited. 

6) The criterion for the growth of the Alfven waves can be derived from 
(16) or (17) 

Vo>~^V A , (20) 

k z 

which means that the velocity of the electron fluid (corresponding to the 
electric current) should be larger than the phase velocity of the Alfven waves 
along the direction of the magnetic field. Otherwise the background plasma 
will attract the free energy of eletromagnetic waves, which is similar to 
Landau damping. In the classical theory for the tearing modes (Furth et al., 
1963) this criterion is the following 

A' = > 0. (21) 

-Dll 

Here e is the width of boundary layer. In fact this criterion means that the 
direction of perturbed electric current should be opposite to the direction of 
the toroidal electric current. 

Note that A' is related to the perturbed parameters, hence it cannot be 
determined by the basic parameters of the plasma. 

7) For the typical parameters of a tokamak plasma, u>, can be estimated 
from Equation (17): if Bo ~ 1 T, n e ~ 10 19 m -3 (a > pe /tt e ~ 1), T e ~ 
iKev, V A ~ 10 6 m/sec, Jo ~ 10 6 A/m 2 , then, w, ~ 10 1 — 10 2 /sec. This 
order of magnitude is comparable with the experiments in tokamaks. In 
the corona the intensity of electric current can not be measured directly. In 
typical solar coronal conditions, n e ~ 10 15 m -3 , Bo ~ 10~ 2 T, T e ~ 10 6 A', so 
V A ~ 10 6 m/sec and cj pe /il e ~ 1 (the similar order as in a typical tokamak 
plasma). The intensity of the electric current in a model coronal current 
sheet is about 10 3 A/m 2 (Martens, 1988), so the value of o>, is from 10° to 
10 1 . The resistivity of the plasma is estimated by the Spitzer’s formula. 

8) Because the Alfven time is an important time scale, the growth of 
Alfven waves may be used to explain many phenomena, such as the quasi- 
periodic oscillations of the solar radio spikes. The author will study this 
subject in a forthcoming paper. 
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Abstract 

The classical theory of cyclotron resonance absorption is discussed in this 
paper. An approximate formula for the cyclotron resonance absorption co- 
efficients is derived from plasma kinetic theory. Finally, the main properties 
of resonance absorption are calculated: the absorption coefficients increase 
rapidly with increasing 6 , the directional angle between the magnetic field 
and the wave-vector, and their values are directly proportional to the elec- 
tron temperature T e and u> pe /Q e . 
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1. Introduction 


The theory of cyclotron resonance absorption is an important part of solar 
radiophysics. Especially solar radio spikes have attracted wide attention 
these years and the electron cyclotron maser instabilities are today’s most 
favoured mechanism for explaining millisecond radio spikes (Benz, 1986; 
Melrose, 1984; Vlahos et al., 1985; Li, 1986; Huang, 1987). But the question 
remains, which modes of the electromagnetic waves can penetrate through 
the corona and will be observed by us? This depends mainly on the magni- 
tude of cyclotron resonance absorption, and there are some different opinions 
on this problem (Melrose, 1984; Sharma et al., 1984; Huang, 1987). 

The classical theory of cyclotron resonance absorption has been estab- 
lished by Stepanov (1958), Gershman (1959, 1960), Zheleznyakov (1973), 
and Ginzburg (1978). The formula that is often used in astrophysics is 
derived by Zheleznyakov (1973), using the single-particle model and Kirch- 
Jio/T’s Law. 

As Ginzburg pointed out (1978), KirchhofTs Law can only be used when 
matter and radiation are in complete thermal equilibrium. When a system 
is unstable (or the waves are amplified), KirchhofTs Law may be used under 
the condition of weak absorption or amplification of the waves, in fact it 
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is only proper for the states in which the emission equals the sum of the 
emissions of all particles. 

But cyclotron resonance absorption is a strong collective plasma effect, 
so the proper treatment is to solve the the kinetic equations (Vlasov’s and 
Maxwell’s equations). Stepanov (1958) and Gershman (1959, 1960) derived 
the kinetic-theory formulae, which can only be used near the regions of res- 
onance. The deviation of the frequency should be smaller than the Doppler 
shift by the thermal motions of the electrons, so this is not absorption under 
the strict resonant condition. In their derivations, the frequency is always 
assumed to be real, and the wave-vector a complex vector, the imaginary 
part of the wave-vector is derived from the plasma dispersion relation and 
is defined as the absorption coefficient. 

In practice the cyclotron resonance absorption is a wave-particle resonant 
interaction in the velocity space, the dimension of which is characterized by 
the Debye Length of the plasma, and it is often much smaller than the 
characteristic length of the plasma uniformity. So the dissipative process of 
electromagnetic waves in a uniform plasma system can only be discussed by 
the imaginary part of the frequency and the group velocity. In a very non- 
uniform plasma, it is difficult to define the group velocity, and the resonance 
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absorption should be studied by solving a boundary value problem for the 
kinetic equations (Ginzburg, 1978). 

On the other hand, when the wave- vector is assumed be a real vector and 
the frequency is complex, for given wave-vector (including the wavelength 
and directional angle), the real and imaginary parts of the frequency will be 
found self-consistently from the dispersion equation. In the opposite case 
(Stepanov, 1958; Gershman, 1959, 1960), the real and imaginary parts of 
the wave- vector should be assumed in the same direction, but this condition 
is not reasonable for an anisotropic plasma. 

The author has derived a formula for cyclotron resonance absorption, 
using the theory of cyclotron maser instability, to explain millisecond radio 
spikes (Huang, 1987). This formula will be simplified further to discuss the 
properties of the cyclotron resonance absorption systematically. In Section 
2 the kinetic-theory formula for absorption coefficients is derived and sim- 
plified. The properties of resonance absorption are discussed in Section 3. 
Section 4 is a summary. 


4 



2. Derivation 


2.1 Physical Analysis 

Cyclotron resonance absorption is a non-collisional absorption mechanism, 
it means a dissipative process of electromagnetic waves under a resonant 
condition. The free energy of the electromagnetic waves is absorbed by 
the ambient electrons, which often have a Maxwellian distribution in the 
velocity space. If the ambient electrons have an anisotropic distribution, it 
will result in the opposite process, the electron cyclotron maser instability. 

Therefore, if the Maxwellian distribution function is used in the formula 
for the growth rates of the electron cyclotron maser instabilities (Wu, 1985), 
instead of the non-thermal distribution of the electrons, the imaginary part 
of the frequency will be negative, which corresponds to the resonance ab- 
sorption coefficient. 

2.2 Derivation 

With the approximations w, C u T and Dfj < D^ (u;,- and u r are the imag- 
inary and real parts of the frequency, Dfj and are the symmetric and 
anti-symmetric components of the plasma dielectric tensor), the exponen- 
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tial decay rates of different modes of the electromagnetic waves in a ‘cold’ 
plasma dispersion relation by cyclotron resonance absorption can be found 
directly from the formula of Wu (1985): 
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On 

Here fc n = a 2 = 221, w pe and fl e are the electron plasma 

and cyclotron frequencies respectively, c is the velocity of light, T e is the 
average electron temperature, m e is the mass of a single electron, k ± and 
are the vertical and horizental wave-vector. In addition, 
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V’r* = 2N 2 (N 2 — e) sin 9 cos 9 , 

V’xy = 25(17 - N 2 sin 2 6), 

V'y* » 2 N 2 g sin 9 cos 9, 

G = 2JV 4 [(1 - C ) 2 £iiil£ + (i_„) cos 2 0]+ 

1 -n 

+4JV 2 [(1 - c ) 2 - 1 - (1 + 0.5 sin 2 9 ■“-]+ 

1-rj 

+2[(1 + i})e 2 -g 2 + 2i} 2 — — — ). 

1 - 1 } 

2 j — 2 

Here s = 1 - , g — ij = 1 - ^r, TV is the refractive index and 

9 is an angle between the propagation vector and the magnetic field. 

3.3 Non- relativistic approximation 

In the corona a typical value of T e is 10 6 A', hence the thermal velocity 
I e ~ 0.01c, and the relativistic factor 7 ~ l. So the resonant condition is 
simplified to 

fUl _ m 

c u r N cos 9' K ’ 

From the intgral formula (Wang et al M 1979) 

f exp(— p 2 t 2 )Jfc(at)J/(6t)t ,+1 dt = 
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Therefore, Equation (2) can be written as 
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2.4 Further Simplification 


• 2 
With the non-relativistic condition, it is easy to prove that p- C 1 (when 

T e = 10 6 Ar, p’ ~ 10 -4 ). If we consider only the lowest power of jj. Equation 

(4) can be simplified to 


£» - ^ 2T e \2 n — 3/2 

fi e “ Guy fi e ' W 2 ' 

(2n — 1 )!! 1 , a 2 .. , , , . 3 

X ~ n! )2 2^ eX P(-p-)((^r + V>yy + 0xy)n + 

4" ""'q “ (ipzz ~~ l})yz 4" Ipxx 4" 0yy 4" V’xy)o^4" 

fc 2 u 2 fc 3 u 3 

+ - V*yz + 0zz)n + (6) 

and T(n + 1/2) = (2n + l)!!7r 1 / 2 /2 n , T(n 4- 1) = n!. 


3. Discussion 


3.1 Directionality 
From Equation (5) we find 


Ui ~ (sin*) 4 ”- 1 . 


so the magnitude of cyclotron resonance absorption is very sensitive to the 
propagation angle 8 of the elctromagnetic waves: the value decreases rapidly 
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with decreasing 6. The same conclusion can also be derived from the classical 
formula (Zheleznyakov, 1973). 

Note that the directionality of the electron cyclotron maser instability 
is quite different for different plasma parameters and different modes of the 
electromagnetic wave, which means that the electromagnetic waves excited 
by the cyclotron maser instability are not always propagating in the direction 
near the vertical to the magnetic field. For example, the directional angle of 
the fundamentals of the fast extraordinary modes is small for typical coronal 
parameters (Huang, 1987). 

On the other hand, the frequencies of the peak values of the cyclotron 
resonance absorption are always in nfi e (n is an integer). But the frequencies 
of the peak values of cyclotron maser instability vary with different modes 
and 6 (Iluang, 1987). 

The calculations of the directionality are summarized in Figures la.b. 


3.2 Effect of T e 


From Equation (5) the absorption coefficients decrease with T t , because 
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which means that the lower the electron temperature, the larger the propa- 
bilitv of the Maxwellian distribution near the origin of the velocity space, 
so the growth rates (depending on an integral along the resonant ellipse) 
decrease. 

In addition, from Equation (3) when U|| = 0 or exp (— 777 1 ) = 1* and at 
the peak values (u; r = nft c ), 


k>i,max 


^ T 2n-3/2 


The results of the relevant calculations are given in Figures 2a, b. 


3.3 Effect of u> pe /Sl e 

It is easy to find from Equation (5) 

So for given cyclotron frequency (or given magnetic field strentgh), in- 
creases with the ambient plasma density, because the more particles are in 
the wave-particle interaction, the stronger the cyclotron resonant absorption 
is (Figures 3a,b). 
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3.4 Absorption of n = 1 


There are some special properties of cyclotron resonant absorption in the 
fast extraordinary modes at n = 1: 

1. When Wpe/n e < 1, the refractive index of the peak value of the 
absorption at n — 1 is found to be in the range of the Slow extraordinary 
modes (N > 1), so the fast extraordinary modes deviate far from the peak 
value of the resonant absorption at n = 1. 

2. Because ~ sin 3 0 at n = 1 from Equation (5), and several terms in 
the brackets of (5) are related to 0, u,' does not vary monotonically with 9. 
The highest peak value is near 6 = 70°, but the efffect of u >,• ~ sin 3 0 is still 
important for the smaller values of 6 (Figure 4) 

3. In Figure 5, the relation of u?,- and T e is similar to that in the Sec- 
tion 3.2, but the peak value of the absorption at n = 1 does not increase 
measurably with increasing T e , because w,-, max ~ Z V 2 - 

4. It is interesting that at higher temperatures, the peak value of absorp- 
tion at n = 1 decreases abruptly to the regions of w; > 0 and a very narrow 
positive peak is found. This means that the high temperature plasma may 
excite the slow extraordinary modes directly (Figure 6). 
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4. Conclusion 


The main point of this paper is that cyclotron resonance absorption is the 
opposite process of the electron cyclotron maser instability, so they can be 
described by the same model (Wu, 1985). This model is more strict than the 
classical single-particle theory, because the single-particle model can not be 
used to explain the strong plasma collective property, such as the resonant 
absorption; and the wave-particle interaction in the velocity space can not 
describe large-scale decay process of the electromagnetic waves directly. Also 
the imaginary and real parts of the wave- vector are not always in the same 
direction. 

A formula for the resonant absorption coefficients is derived from plasma 
kinetic theory under the approximation of u;, < w> and D,“ < and the 
main properties of the cyclotron resonant absorption are: 1. w, increases 
rapidly with increasing 6, 2. the peak value of u>,- increases rapidly with 
increasing T e , 3. w,- is directly proportional to the parameter u,* pe /n e . 4. 
There are some special properties for the absorption of the extraordinary 
modes at n = 1, which will be studied further by the author. The formula is 
also valid for the absorption of the ordinary modes and Cerenkov absorption 
at n — 0. 
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Note that cyclotron resonance absorption will result in anisotropy of a 
Maxwellian distribution, which is the basic difference with collisional ab- 
sorption. So it is a possible mechanism for the secondary acceleration, and 
to exite a new electromagnetic emmision. The average kinetic energy and 
density of the accelerated particles can be estimated by the theory in this 
paper. This subject will be discussed in a forthcoming paper. 
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Figure Captions 

Fig. la.b: The relation of the second and third absorptions of the fast 
extraordinary modes with the angle 0 (B = 200 Gauss , n e = 10 s c7n~ 3 , 
v; = Ufk). 

Fig. 2a,b: The relation of the second and third absorptions of the ex- 
t inordinary modes with T e (B = 200 Gauss, n e = 10 8 cm -3 , 0 = 60°). 

Fig. 3a,b: The relation of the second and third absorptions of the ex- 
traordinary modes with u> pe /fl e (B = 200 Gauss, 6 = 60°, T e = 10 6 A'). 

Fig. 4: The relation of the absorptions of the extraordinary modes at 
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» = 1 with 6 (B = 200 Gauss , n e = 10 8 cm* 3 , T e = 10 6 A'). 

Fig. 5: The relation of the absorptions of the extraordinary modes at 
n = 1 with T e ( B = 200 Gauss, 6 = 80°, n e = 10 8 cm" 3 ). 

Fig. 6: The relation of the absorptions of the extraordinary modes at 
n = 1 with Wpe/fie ( B = 200 Gauss, 6 = 80°, T« = 10 6 A"). 
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ABSTRACT 

In this paper an extension of the classical theory of the tearing mode Instabilities Is dis- 
cussed. It Is suggested that the periodicity of the tearing modes should be considered, and 
that the induced electric field will be available to accelerate the energetic particles to 
excite the millisecond radio spikes In solar flares. Moreover, the period of the tearing 
modes agrees with the timescale of a single spike, and the characteristic time of the exponen- 
tial growth of the tearing modes is comparable with the timescale of a group of spikes. So 
the evolution of the millisecond radio spikes may be explained by the quasl-perlodlclty of th« 
accelerating electric field or the. energetic particles. 


INTRODUCTION 

The particle acceleration mechanism Is an Important subject of solar flare physics; It deals 
with the process of magnetic energy conversion into the kinetic energy of particles, then to 
excite various electromagnetic radiations. Including hard X-rays, Type III bursts, mllllseconc 
radio spikes, microwave radiation, etc. 

It is noted that the source of radio spikes is close to the primary acceleration region; that 
Is to say, the radio spike Is a direct signature of the energetic particles /I /. On the other 
hand, the research of acceleration mechanism may deepen the understanding of emission mecha- 
nism of the radio spikes. As de Jager /2/ pointed out In his recent review, there are many 
kinds of acceleration mechanisms, but for an impulsive process (such as the radio spikes), th« 
acceleration by electric field Is the only possible mechanism. Solar radio spikes have 
attracted wide attention these years, and the electron cyclotron maser Instabilities are 
today's most favored mechanism for the millisecond radio spikes in declmetrlc range /!/, many 
characters of which are explained by this theory: high bright temperature, narrow bandwidth, 

directionality, polarisation, harmonic structure, etc. /3, 4,5, 6, 7, 8/. 

• 

However, many problems are still unclear. Including the evolution of spike emission or the tit 
structure of a single spike and a group of spikes. This problem may be related to the proper- 
ties of energetic particles or acceleration electric field /9,I0,II/« So the question Is to 
find a particle acceleration mechanism, which could accelerate the energetic particles needed 
to excite the millisecond radio spikes and explain the evolution of radio spikes. 

Tearing Mode Instabilities and their Induced Electric Field 

Restriction of the classical theory. The resistive tearing mode instability Is one of the 

most Important large-scale MHD instabilities. It may result In the disruption of plasma in 
tokamaks, so the tearing mode Instability Is studied Intensively in laboratories and in 
theories of plasma physics. 

In plasma astrophysics, Heyvaerts /12/ and Spicer /13/ also discussed the acceleration 



mechanism of tearing modes. It is known that when resistivity Is considered In the MHD 
equations, the hypothesis of frozen-in field line is no longer correct, and then the re- 
sistive MHD instabilities may be excited by the electric current (the free energy of magnetic 
field). Because of the shear of the magnetic field and the periodic boundary conditions, 
these instabilities can only be excited near some spatial magnetic surfaces, which are often 
called "rational surfaces". Therefore, in the classical theory of tearing mode instabilities 
/14/, a WKB approximation is used: the resistivity is considered only in a boundary layer, 

and the ideal MHD equations are still valid in other regions. Then the growth rates of the 
tearing modes (the real parts of eigenvalues) are derived from the continuity boundary con- 
dition. 

But in the classical theory, the imaginary parts of eigenvalues are not taken into consider- 
ation. In fact the imaginary parts of eigenvalues describe the periodic oscillation of 
tearing modes, and the real parts of eigenvalues mean an exponential growth of the amplitude 
of oscillations. The author has shown theoretically /I5/ that the evolution of nonlinear 
tearing modes is a quasi-periodic process. So its Induced electric field and accelerated 
particles must vary with the same period, which may be used to interpret the evolution of 
radio spikes. 

Derivations . The resistive MHD equations are / 16/ s 


3Y/3t + (u # V a )Y - nVx Y - E 
3U/3t + Ou^JU - x Vi(7i f)l 

u - V x $ x z 
B - -rV x Y x z ♦ z 
J 
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Vx Y 


U - -V± * 


( 1 ) 

( 2 ) 

( 3 ) 

(4) 

( 5 ) 

( 6 ) 


Here Y and ♦ are the potential functions of magnetic field B and velocity of magnetic fluid 
u respectively, n is the resistivity of the plasma, E is the electric field related_to the 
electric current along the magnetic field lines; the parameter S • T p/ T u* T g " • 
t h * Wy^p^) 1 2 /B, r - a/R, R and a mean the large and small radius of a toroidal plasma, 

fj is the average resistivity of the plasma, p n is the mass density of the magnetic fluid, 
p is the magnetic conductivity in vacuum. Equations (l)-(6) are dimensionless: L/a, 

t?t R , n/n. B/B s . J x /(B r /P o R), 7/« J B z . 

Then it is postulated that the growth of tearing modes depends mainly on the properties near 
the corresponding rational surface, and one is mainly interested in the evolution of tearing 
modes; hence, the solutions of ( 1 ) — ( 2 ) will be averaged in the boundary layer. 


In the plane slab approximation, x is the radial coordinate, and y is the pololdal coordinate. 
Let : 
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Then the linear approximation of equations (l)-(2) can be written: 
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Her. < > means an average in the boundary layer, and 

of t and x respectively. Let: 


express the derivatives 


X - - toWH/N 
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( 16 ) 


Y - V (17) 

T “ 0) t (18) 
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U) 2 o _ -CN/AH (1*) 

Equations (9)-(10) can be written as: 

- X + Y - 0 (20) 

Y + X + Y - 0 (21) 

8 - D/0) A < 0 (22) 

o 

It is easy to find the solutions of equations (20)-(2I): 

X - c^exp(-gT/2) cost (23) 

Y - c 2 exp(-gx/2) sint (24) 


c 1 and are the integral constants, and the condition |g| < 1 is used. 

Discussion 

Comparison with the classical theory . From equations (23), (19), (11)-(15) and (3)-(6), the 
growth rates of the tearing modes are: 

- ga) o /2 - kn < J zi >/< B xl > (25) 

If the parameter A' and the width of boundary layer h are defined: 

A* - <B \ *> h/<B ,> (26) 

x 1 xl 

Then from Ampere's law: 

< 3 zl> * A'<B xl >/kh (27) 

i the growth rate in (25) may be changed into: 


- gu) o /2 - nA'/h (28) 

This is the result of Bateman /I7/ v and because: 

h -v S -2/5 (A') 1/5 (29) 

- «» 0 /2 ' V 3 ' 5 V 2 ' 5 W )‘" ( 30 ) 

This is consistent with the classical theory /14/. On the other hand, the imaginary part of 
eigenvalue can be derived from (19), (11)-(15), and (5)-(8): 


In fact, *(31) is the dispersion equation of tearing modes. 

Quantitative analysis . In the typical coronal condition: R^l0 4 km, a'w i0 2 km. r - a/R 
* 10 it can be estimated from the MHD equations /16/t B %r, u .%r, J ^r , B ,%r , 
J 2 ' 0 '^r, U x - \\ x u| ^ r ..., the period of tearing mod£2 equals: * 

T - t R r/s Vf^r % 10‘ 2 sec. 


This is comparable with the typical duration of a single spike. Moreover, the characteristic 
time of the exponential growth of the tearing modes is: 

T R /gu> 0 + T r ^ 1 sec (33) 

This is of the same order as the typical duration of a group of spikes. Otherwise, the phase 
velocity of tearing modes may be estimated by (31): * 

u o 2 /k 2 ^ 2 S 2 / (T r 2 r 2 ) ^ R 2 / tjj 2 - u A 2 ( 3 *) 


This Is Just the value of Alfven velocity 


Note that the parameter |g| % r « 1. 

Properties of the Induced electric field . The typical order of magnetic field B in the 
corona is io 2 G - l<f 2 T; therefore B \ r 2 B \ 10~ 4 T, and u a 10 6 m/sec /18/. Then the in- 
duced electric field can be derived By Faraaay’s lav: 

E. ^ uj B./k ^ 1 v/m. 
l o l 

This electric field can accelerate an electron to tens of KeV during a period of the tearing 
Bodes and satisfy the requirement of the radio spikes for the energetic particles. 

Spicer /13/ est lasted the Induced electric field by the exponential growth rates: 

E. 't* gu) B/k % 10 2 v/b. 
i o 

This value is two orders of magnitude gmaller than that in (35). On the other hand, the 
characteristic length of E^ Is: L'Va'VlO -km, which equals the distance passed by an electron 
at the velocity of light during a Billisecond. The energy gain of an electron over this 
length scale is Indeed several tens of KeV. 

INTERPRETATION OF THE TIME STRUCTURE OF THE RADIO SPIKES 

Pros the Boaentum equation of an electron: 

a^dv/dt - eE^ cos w o t (36) 

It is easy to get: 

v - eE./m ci) sin CD t + "v . (37) 

i e o o o 

o 

Here the aaplltude of velocity 0) is about equal to 10 m/sec, which Is of the same 

order of magnitude as the velocity of Sight in vacuum. Moreover, in a given direction the 
energetic electron beam is not uniform but of a periodic intermittence. 

It Is noticed that the growth of electron cyclotron Baser instabilities depends on the 
energetic electron beam, so the intermittence of the beam will make the spike emissions decay 
rapidly to the background level. Otherwise in (31) the frequency u) Is proportional to B . 
That Beans the acceleration time decreases with Increasing the magnetic field, which may 20 
be used to explain why the duration of radio spikes decreases with increasing the frequency 
/l/ For example, the duration of radio spikes In the Betric range Is one or two orders of 
magnitude larger than that in the declmetrlc range. 

SUMMARY 

In this paper an extension of the classical theory of tearing mode instabilities is discussed 
It is suggested that the periodicity of tearing modes should be considered, and that the 
Induced electric field will be available to accelerate the energetic particles that excite 
the radio spike emission. Moreover, the period of the tearing modes corresponds with the 
tiae scale of a single spike, and the characteristic time of the exponential growth of the 
tearing modes is comparable with the time scale of a group of spikes. Therefore the 
acceleration mechanism of the periodic tearing modes gives a reasonable model to describe the 
Billisecond radio spikes in solar flares. 
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